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ABSTRACT 


This  thesis  is  a  theoretical  study  of  conventional  and  unconven¬ 
tional  potential  models  for  proton  inelastic  scattering.  The  unconven¬ 
tional  potential  results  from  a  reduction  of  the  Dirac  equation,  with 
scalar  and  vector  complex  nuclear  potentials,  to  a  Schrodinger-type 
equation.  At  intermediate  energies  the  potential  is  found  to  have  an 
unconventional  shape,  namely  the  Wine  Bottle  Shape.  This  potential  is 
generalized  to  an  extended  optical  potential  appropriate  for  a  treat¬ 
ment  of  inelastic  excitations  of  collective  states.  The  generalized 
potential  is  found  to  have  a  deformed  spin  dependent  interaction  of 
the  full  Thomas  form.  Comparisons  are  made  between  predictions  for 
inelastic  scattering  observables  calculated  using  this  potential  and 
those  based  on  a  non-relativistic  Schrodinger  equation  based  optical 
potential  of  the  Standard  form.  Results  of  these  comparisons  for  both 
cross-section  and  analyzing  power  will  be  discussed. 

The  effects  of  using  the  deformed  spin-orbit  interactions,  the 
deformed  Coulomb  interaction,  relativistic  versus  non-relativistic  kine¬ 
matics  and  equal  deformation  parameters  (3)  versus  equal  deformation 
lengths  (3R)  are  investigated. 

The  main  conclusion  of  this  work  is  that  the  unconventional  poten¬ 
tial  (Dirac  model)  predictions  for  proton  inelastic  scattering  cross- 
section  and  polarization,  for  some  low-lying  collective  states,  are 
either  comparable  or  superior  to  those  of  the  conventional  potential 
(Standard  model) .  We  have  found  that  the  deformation  length  ftR,  deter¬ 
mined  from  a  comparison  of  experimental  and  calculated  cross-sections, 
decreases  with  increasing  projectile  energy. 
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I.  INTRODUCTION 


One  of  the  most  important  concerns  in  nuclear  physics  is  to  under¬ 
stand  the  structure  of  the  nucleus  as  well  as  the  mechanisms  of  nuclear 
reactions.  In  practice  one  can  find  the  properties  of  nuclei  by  appro¬ 
priate  scattering  and  reaction  measurements.  However,  the  interaction 
of  nuclear  projectiles  with  nuclei  is  a  complicated  many-body  problem. 
Because  of  this  complexity  one  often  relies  on  simplified  models  to 
describe  the  phenomena. 

One  of  these  phenomenological  approaches  is  the  optical  model 
potential  based  on  Schrodinger  equation  with  Woods-Saxon  form  factors, 
which  seeks  to  describe  the  interaction  of,  say,  a  nucleon  with  a  tar¬ 
get  nucleus  by  means  of  a  complex  potential.  This  model  has  been  suc¬ 
cessful  in  describing  proton  elastic  scattering  at  low  energies  and 
few  hundred  MeV  energies.  We  will  call  this  approach  "The  Standard 
Model"  throughout  this  work. 

In  recent  years,  another  optical  model  potential  based  on  the 
Dirac  equation  has  been  constructed  from  a  Lorentz  scalar  term  and  a 
time-like  Lorentz  vector  term.  It  has  been  found  that  this  model  is 
successful  in  describing  the  proton  elastic  scattering  at  intermediate 
energies  (Ar81,  Ar82 ,  Co81) .  One  of  the  interesting  results  of  this 
model  is  the  discovery  of  the  "wine  bottle  shaped"  real  central  poten¬ 
tial.  It  is  found  that  when  one  reduces  the  Dirac  potential  to  "Schro¬ 
dinger  equivalent  potentials"  that  the  real  central  potential  turns 
repulsive  in  the  interior  at  energies  between  100-200  MeV  with  attrac¬ 
tion  in  the  outer  region  (Ja79,  Ja80) .  It  becomes  progressively  more 
repulsive  as  the  energy  increases.  This  is  in  contrast  with  the  standard 
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Woods-Saxon  potentials  where  the  real  central  potential  weakens  but  re¬ 
mains  totally  attractive  as  the  energy  increases,  then  turn  repulsive 
at  some  energy.  (In  this  work,  we  will  call  the  latter  optical  model 
approach  MThe  Dirac  Model" .) 

The  parameters  of  both  potentials  can  be  treated  completely 
phenomenologically  and  thus  determined  from  the  data,  or  the  number  of 
free  parameters  can  be  reduced  by  invoking  nuclear  structure  informa¬ 
tion  on  charge  and  matter  density. 

Both  potentials  have  been  generalized  to  describe  the  inelastic 
scattering  processes  leading  to  the  excitation  of  collective  surface 
modes  of  target  nuclei.  In  this  generalization  the  normal  spherical 
potential  is  extended  by  including  terms  that  explicitly  contain  dyna¬ 
mical  nuclear  coordinates  describing  the  displacement  of  the  nuclear 
surface.  This  is  usually  achieved  by  letting,  R+R+a(r),  where  R  is 
the  radius  parameter  of  a  typical  spherical  optical  potential  a(r)  is 
the  displacement  of  the  nuclear  surface  in  the  direction  r.  All  parts 
of  the  potential  are  then  deformed  and  expanded  in  powers  of  a(r)  and 
only  terms  to  first  order  in  a(r)  are  retained  (B164) .  The  calcula¬ 
tions  of  the  inelastic  scattering  cross-section  and  polarization  are 
done  in  the  framework  of  the  distorted  wave  Born  approximation. 

In  a  recent  work  Satchler  (Sa83)  has  studied  the  inelastic  cross- 
section  for  low-lying  collective  states  in  ^Ca  anci  for  proton 

energies  near  200  MeV.  His  study  was  based  on  extending  two  spin  inde¬ 
pendent  optical  potentials,  that  provide  the  same  elastic  scattering 
cross-section,  to  describe  the  inelastic  scattering  section.  One  of 
these  potentials  is  a  potential  of  standard  Woods-Saxon  shape,  while 
the  other  is  of  a  wine  bottle  bottom  shape  generated  through  the  addition 
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of  a  derivative  term  to  the  potential.  Satchler  did  not  include  spin- 
orbit  interactions  in  his  analysis  of  inelastic  scattering. 

In  the  present  work  we  will  carry  out  a  similar  investigation  of 
the  consequences  of  the  unconventional  shape  of  the  nuclear  optical  po¬ 
tential  at  intermediate  energies  for  proton  inelastic  scattering.  Our 
approach  is  based  on  comparing  the  calculations  resulting  from  both  the 
standard  and  the  Dirac  models.  Our  investigation  differs  from  that  of 
Satchler  in  one  additional  respect,  namely,  that  the  spin-orbit  effects 
are  fully  taken  into  account. 

In  chapter  II,  we  describe  the  standard  optical  potential  model 
for  proton  elastic  scattering.  In  chapter  III,  we  describe  the  extended 
standard  optical  model  for  proton  inelastic  scattering.  In  chapter  IV, 
we  describe  the  Dirac  optical  model  for  proton  elastic  scattering.  In 
chapter  V,  we  describe  the  extended  Dirac  optical  model  for  proton  in¬ 
elastic  scattering.  In  chapter  VI,  fits  to  elastic  scattering  data  from 
both  models  and  comparisons  of  the  conventional  and  the  unconven¬ 
tional  potentials  will  be  shown.  In  chapter  VII,  theoretical  predic¬ 
tions  of  the  present  work  for  inelastic  scattering  are  compared  to 
experimental  data,  and  to  Satchler Ts  predictions  when  possible.  A 
conclusion  is  given  in  chapter  VIII. 


II.  THE  OPTICAL  POTENTIAL  MODEL 


If  a  nucleon  enters  a  nuclear  field  of  a  target  nucleus,  it  will 
interact  with  the  target  in  its  path  and  will  be  either  absorbed  or 
emerge  on  the  other  side  after  some  considerable  buffeting.  At  each 
encounter  it  may  change  its  direction  and  spin  state,  and 
several  such  collisions  may  take  place  as  the  incident  particle  crosses 
the  nucleus.  The  interaction  of  a  single  nucleon  with  a  nucleus  is  a 
complicated  many -body  problem.  The  optical  model  replaces  the  many- 
body  interaction  by  a  complex  nucleon-nucleus  potential.  In  optical 
model  phenomenology  the  potential  is  written  in  terms  of  certain  model 
parameters  which  are  determined  from  comparison  with  experimental  data. 

A  formal  theory  is  needed  to  relate  these  model  parameters  to  the  pro¬ 
perties  of  the  many-body  system. 

II. 1  Formal  Theory  of  the  Optical  Model  Potential: 

We  will  summarise  here  the  formal  theory  of  the  optical  model  deve¬ 
loped  by  Feshbach  (Fe58)  as  it  was  described  by  Jackson  (Ja70) . 

If  we  consider  a  nucleon  interacting  with  a  nucleus,  the  total 
Hamiltonian  of  the  system  could  be  written  as 

H  =  Hq  +  H(£)  +  V(?,0  (II. 1) 

where  Hq  is  the  kinetic  energy  operator  of  the  relative  motion  of  pro¬ 
jectile  and  target,  H(£)  is  the  internal  Hamiltonian  of  the  target  nu¬ 
cleus  which  has  co-ordinates  £,  and  V  is  the  interaction  potential  bet¬ 
ween  the  projectile  and  the  target  nucleons.  The  Schrodinger  equation 
we  have  to  solve  here  is 

HT(?,£)  =  ET(?,0  (II. 2) 
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where  ¥(r,0  is  the  total  wavefunction  of  the  projectile-target  system. 

The  target  states  (Q  are  solutions  of  the  equation 

A 

H(0«x(C)  =  ^(O  (II. 3) 


where  e  is  the  energy  eigen  value  of  the  target  Hamiltonian  H(£)  which 


is  assumed  to  be  hermetian.  The  states  are  completely  antisymme¬ 

tric  and  form  a  complete  orthogonal  set  with  the  ground  state  denoted  by 
$  (£) .  The  total  wavefunction  T  may  be  written  as  an  expansion  in 


=  l  ^(?)^(?) 


X 


X 

(II. 4) 


where  i|j^(r)  describes  the*  motion  of  the  projectile  relative  to  the  target . 

The  elastic  channel  in  the  above  expansion  is  represented  by 
^o(r)$o(0 .  If  we  define  P  as  an  elastic  channel  projection  operator: 

p¥(?,£)  =  i|>o(?)*o(£).  (II.  5) 

and  if  we  introduce  another  projection  operator  Q  such  that  it  projects 
out  all  states  except  the  elastic  channel,  we  have 


PQ  =  QP  =  0 

then  one  can  write  the  following 
Q¥(r,£)  =  (1-P)Y(?,£) 

=  -  Pl-Cr.O 

=  -  ipo(?)$oce) 


(II. 6) 


(II. 7) 


We  can  construct  these  projection  operators  as  follows 


. 
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p  =  |$o(C)x$o(5)|  (ii. 8) 

and 

Q  =  l  IV<*,|  •  (II-9) 

A#>  A 

Equation  (II. 2)  can  now  be  rewritten  as 

(E  -H)(P+Q)T  =  0  (II. 10) 

if  we  multiply  on  the  left  by  P,  we  get 

(E-PHP)P'i'  =  (PHQ)Q'i'  (II. 11) 

on  the  other  hand*  multiplication  by  Q  yields 

(E-QHQ)QT  =  (QHP)PT  (11.12) 

where  we  have  used  the  properties  of  the  projection  operators  as  given 
by  equation  (II. 6).  We  can  rewrite  equation  (11.12)  as  follows 

=  rdW QHP  **  •  (II-13) 

If  there  are  open  non-elastic  channels  it  is  necessary  to  replace  E-QHQ 
in  equation  (11.13)  by  E-QHQ+ie,  in  order  to  ensure  that  QT  has  only 
outgoing  waves  in  these  channels.  If  we  substitute  equation  (11.13) 
into  equation  (II. 11)  we  get  the  following  equation  for  P¥ 

(E  -  PHP  -  PHQ  9 — 9—  QHP)Pf  =  0  .  (11.14) 

.fcj  —  QriQ 

Since  we  are  interested  only  in  elastic  scattering  we  have  to  reduce 
equation  (11.14)  to  an  equation  for  ^(r).  This  can  be  done  by  multi¬ 
plying  from  the  left  by  <$Q|  and  integrating  over  the  target  co-ordinate 
£.  Then  one  can  write  equation  (11.14)  as  follows 

<$o|E|PT>  -  «I>o  I  PHP  I  PT>  -<$JPHQ  e  -  QHP  |  PT>  =0  (11.15) 


and  hence 


. 


7 


EifJ  (r)  -  H  ifj  (r)  -  £  i p  (r)  -  <$  |v(r,£)|$  (r) 

o  o  o  o  o  o'  o  o 


lPHQ  eTT^htT  QHP  |  *>«„(?)  =  0 


“  <$0 I PHQ  E  _  qhq 


(11.16) 


the  last  term  can  be  reduced  as  follows: 


<$  I PHQ 
o ' 


l  <*  |  PH  |  «-><$.  | 

W  °  X  A 

I  <$ Jh  +  H(C)  +  V(r,0  |$  ><$.  | 
A#>  °  0  A  A 

l  <$  |v(r,£)  | $.><&.  | 

A#0  0  X  A 

<4>o|v(?,5)Q 


(11.17) 


where  the  first  two  terms  drop  out  because  of  orthogonality,  since 


<$  lH  |  $,  >  =  H  <$  |$,>  and  <$  |h(£)|<1\>  =  E,<$  I  $>  >  (11.18) 

o1  o'  A  o  o'  A  o1  -  1  A  A  o'  X 

if  we  assume  that  £q,  the  energy  of  the  ground  state,  is  zero,  and  if 

we  substitute  equation  (11.17)  into  equation  (IP. 16) ,  we  can  write  the 

latter  as  follows 

[E  -  Hq  -  <$q | V | 0q>  -  <$o  j VQ  Qv|$o>]^(b  =  0  (11.19) 

so  we  can  formally  define  a  generalized  optical  potential  as  follows 

U(r)  =  <$q | V | $q>  +  <$o|vQ  — ^  QV|$0>  (11.20) 

then  we  can  write  equation  (11.19)  as 

(E  -  Hq  -  U (r) )^q (r)  =  0  .  (11.21) 

The  leading  term  in  equation  (11.20),  <$o | V | $q> ,  yields  the  real  parts 
of  the  potential,  and  hence  one  can  approximate  U  as 


(11.22) 


where 
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V  =  V(r,£)  =  l  v. (r  - r . )  (11.23) 

i=l  1  1 

then  one  can  write  equation  (11.22)  as 

u  =  <$  iy  v. (?-?.)i*  > 

o  1  T  1  1  1  O 

1 

=  [dr'  v(r-r')  5^ -r ')  |  $o>  .  (11.24) 

J  i 

But  by  definition 

p (P)  =  «S>  IT  6 (r .  -rf)|$  >  (11.25) 

i 

where  p  is  the  matter  distribution  density,  and  for  short  range  v  one  may 
use  the  6-interaction  and  get  the  following 

U  =  const  p(r).  (11.26) 


Thus  the  potential  follows  the  density  in  this  approximation,  and  hence 
one  is  justified  in  using  the  Woods-Saxon  functions  at  least  at  low 
energies. 

This  model  could  be  extended  to  describe  the  inelastic  scattering. 
This  is  usually  done  by  deforming  the  optical  potential.  Since  the 
potential  follows  the  density,  one  can  take  the  spherical  density  and 
then  deform  the  surface.  This  will  be  discussed  in  Chapter  III. 


II. 2  The  Phenomenological  Optical  Model 

This  approach  makes  the  optical  model  one  of  the  simplest  and 
most  successful  of  nuclear  models.  The  phenomenological  optical  potential 
which  is  used  usually  to  describe  nucleon -nucleus  elastic  scattering 
has  the  following  form 


■ 


' 
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U(r)  =  Vc(r)  -  VqFo  -  i(W  -  ^  -)F. 


m  c 

If  . 


where 


Fx  ^  [1  +  exp 


[V  — —  F  +iW  —  —  F  ]G*L 
L  SO  r  dr  SO  SOI  r  dr  SOIJ 


r  -  R  ^  -1 
x 


(11.27) 


(11.28) 


1/3 

is  the  Woods-Saxon  function  and  =  r^A  ,  in  general  the  radius 

and  the  diffuseness  parameters  r  and  a  are  different  for  different  parts 

x  x 

of  the  potential.  V^(r)  denotes  the  Coulomb  potential  due  to  a  uni¬ 
formly  charged  sphere  of  radius  R^  and  charge  Ze,  which  is  necessary 
for  proton  scattering  to  take  account  of  the  interaction  of  the  incident 
proton  with  the  charge  distribution  due  to  the  protons  in  the  nucleus. 

It  has  the  form 

2 


=  ft1  C3  -  (^>2) 

c  c 


r  <  R 


Ze 

r 


r  >  R 
—  c 


(11.29) 


The  quantity  Vq  denotes  the  depth  of  the  real  central  potential,  W 
denotes  that  of  the  volume  part  of  the  central  imaginary  potential, 

denotes  that  of  the  surface  peaked  imaginary  potential,  and  the  spin- 
orbit  term  is  of  the  Thomas  form  with  V  as  the  strength  of  its  real 

bU 

part,  and  W  that  of  the  imaginary  part.  The  spin-orbit  potential  is 
necessary  for  particles  with  non-zero  spin  because  it  is  needed  to  ex¬ 
plain  the  elastic  polarization  of  the  scattered  particles,  and  also  it 
was  found  that  it  improves  the  calculated  cross-section  at  backward 
angles  (Sh68). 

A  phenomenological  optical  potential  with  suitably  adjusted 


' 


■ 
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parameters,  is  able  to  give  a  good  overall  account  of  the  experimental 
data  on  elastic  scattering.  There  are  thirteen  different  parameters  in 
the  potential,  they  could  be  varied  in  search  for  a  fit  to  the  experi¬ 
mental  data.  In  practice  one  usually  is  restricted  to  a  lesser  number 
of  parameters  to  be  in  search,  this  depends  on  the  proton  incident 
energy.  At  low  energies,  W  is  usually  zero,  and  at  intermediate 

bU 

energies,  is  zero,  so  in  each  case  we  exclude  three  parameters  from 
the  search.  These  optical  model  parameters  are  generally  smooth  and 
slowly  varying  functions  of  energy  and  atomic  weight  of  the  target 
(Na81). 


The  final  parameters  are  determined  usually  by  minimum  values  of 
the  total  Chi-square,  and  sometimes  by  putting  a  constraint  on  the 
search  to  give  a  reasonable  reaction  cross-section.  The  total  Chi- 
square  is  given  by 

X2  =  xl  +  x2  (11.30) 

a  p 

where  a  and  p,  denote  cross-section  and  polarization  respectively,  and 


Na 

l 

i=l 


aex(6i>-°th<9i> 


Act  (9i) 
ex 


(11.31) 


where  aex(9i)  >  O^Oi)  and  Acr^CQi)  are  respectively  the  experimental 

theoretical  and  experimental  error  values  for  the  cross-section  at  a 

2  2 

center  of  mass  angle  0i.  A  similar  expression  holds  for  Xp»  the  X 
associated  with  the  polarization  data. 

The  optical  model  has  been  useful  in  interpreting  nuclear  scatter¬ 
ing  at  low  and  intermediate  energies.  In  particular,  the  occurrence  of 
giant  resonances,  the  total  elastic  scattering  cross-section,  the  total 
reaction  cross-section,  polarization,  and  the  angular  distribution  of 


■ 

■ 
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3 

elastically  scattered  protons,  neutrons,  deuterons,  a-particles,  He, 

3 

H  and  heavy  ions,  can  be  understood  in  terms  of  the  optical  model. 

The  model  is  useful  also  in  describing  the  inelastic  scattering  and 
various  rearrangement  collisions  involving  transfer  of  one,  two  or 
more  particles.  In  general  the  model  is  used  extensively  in  studying 
nucleon-nucleus  scattering,  light  ions  scattering,  heavy  ions 
scattering  and  other  scattering  processes.  Because  of  these  facts  one 
can  say  that  the  optical  model  has  produced  a  considerable  amount  of 
information  about  nuclear  structure  (Ho63) . 

We  have  used  the  optical  model  in  describing  the  proton  nucleus 
elastic  scattering  from  different  target  nuclei  at  various  energies. 
The  results  will  be  analyzed  and  discussed  in  chapter  VI. 


III.  THE  NON-RELATIVISTIC  EXTENDED  OPTICAL 


MODEL  "THE  STANDARD  MODEL" 


If  a  nucleus  is  bombarded  with  a  proton  many  things  can  happen. 

The  most  common  occurrences  are  the  following:  the  proton  can  elasti¬ 
cally  scatter  from  the  nucleus  leaving  the  latter  in  its  ground  state; 
the  proton  can  give  up  some  of  its  kinetic  energy  leaving  the  nucleus 
in  a  bound,  excited  state;  or  else  the  proton  can  break  the  nucleus  if 
it  has  high  kinetic  energy.  In  this  chapter  we  will  discuss  the  direct 
proton  nucleus  inelastic  scattering  where  the  proton  leaves  the  nucleus 
in  an  excited  state.  We  will  consider  only  the  excitations  of  the  low 
lying  collective  states  namely  the  vibrational  and  rotational  states. 
The  discussion  will  be  represented  in  the  framework  of  the  distorted 
wave  Born  approximation. 

Consider  the  total  Hamiltonian  for  a  nucleon  interacting  with  a 
target 


H  =  Hq  +  H(0  +  V(?,D 


(HI.l) 


where  Hq  is  the  K.E.  of  relative  motion,  H(£)  is  the  target  Hamiltonian, 
and  V  is  the  interaction  potential  of  the  incident  nucleon  with  the 
target  and  it  could  be  written  in  the  following  fashion 


A 

V  =  Y  v 


(HI. 2) 


where  A  is  the  target  mass  number,  and  v_.  is  the  interaction  of  the 
projectile  with  the  jth  nucleon  in  the  target. 

If  we  take  T  as  the  total  wavef unction  of  the  target  projectile 
system  such  that 


(H  -  E)T  =  0 


(HI. 3) 
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Then  the  exact  transition  matrix  for  an  inelastic  process  from  an  ini¬ 
tial  state  i  to  a  final  state  f  is  given  by 

Tfl  =  «|>(kf)$f  |v|^+)>  (III. 4) 

where  V  is  the  interaction  potential  in  the  outoing  channel,  the  sub¬ 
script  (+)  refers  to  outgoing  boundary  conditions,  is  the  wavef unc¬ 
tion  of  the  target  in  the  final  state,  and 
-*-kf  * r 

4>f(kf)  =  e  n  (III. 5) 

•> 

where  T]  is  the  nucleon  spin  function,  and  is  the  wave  vector  of  rela¬ 
tive  motion  in  the  final  state  (Ja70) .  Let  us  introduce  a  spherical 
optical  potential  in  an  artificial  way: 

T_.  -  <cj) (k  )$  |  (V-U  )  +U  .  (III. 6) 

ii  t  t  o  o  i 

The  optical  model  distorted  waves  are  defined  by 

(H  +  U0  “  E)X  =  °*  (III. 7) 

Using  the  Gell-mann-Goldberger  relation  for  scattering  from  two 
potentials,  we  get  the  exact  T-matrix 

T. .  -<xi-)*Jv-U  |^+)>  +  <xi-)$Ju  |*.<f>(k.)>  .  (III. 8) 

ri  rr  oi  troll 

Because  Uq  is  spherical 

«£>f |U  |$  >  =0  if  i  ^  f  .  (III. 9) 

Then  for  inelastic  scattering 


(III. 10) 
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which  is  still  exact.  We  do  not  know  ,  so  we  have  to  use  approxi¬ 

mate  methods,  for  example,  the  distorted  wave  Born  approximation  (DWBA) . 
It  is  assumed  that  the  dominant  part  of  ¥  represents  elastic  scattering. 
This  is  supported  by  the  observation  that  in  most  cases,  elastic  scat¬ 
tering  is  the  dominant  process.  If  Uq  is  the  optical  potential  that 
describes  the  elastic  scattering  of  the  projectile  on  the  target  then 
we  may  write 

.  (III. 11) 

1  i  i 

Where  this  is  expected  to  be  a  good  approximation  in  the  asymptotic  re¬ 
gion  but  may  be  bad  in  the  interior  (Ja70) .  The  right  hand  side  of 
equation  (III. 11)  is  also  the  leading  term  in  a  Bom  series  expansion 

T  .  =  |<D.X.(+)>  .  (III. 12) 

fi  f  f 1  o 1  1  1 

We  would  like  to  apply  this  expression  to  the  inelastic  scattering 
leading  to  the  excitation  of  collective  states.  The  collective  model 
of  the  nucleus  allows  for  a  very  simple  excitation  of  the  nucleus  through 
a  change  in  the  collective  state  of  the  system.  The  states  most  strong¬ 
ly  excited  by  inelastic  scattering  are  those  recognized  as  involving 
collective  motion  of  some  kind,  vibrational  or  rotational.  They  usual¬ 
ly  have  the  highest  cross-sections,  this  is  because  of  the  strong  over¬ 
lap  between  the  initial  and  final  states  of  the  nucleus.  The  excita¬ 
tion  of  these  states  is  most  commonly  represented  by  means  of  a  macro¬ 
scopic  model  (the  extended  optical  model) .  The  first  matrix  element  on 
the  right  hand  side  of  equation  (III. 12)  may  be  approximated  as 

<$f  =  <4>f  |UD(r)  (III. 13) 

-y 

where  U^(r)  is  the  deformed  optical  potential.  This  can  also  be  written 
as 
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M?)  =  U  (r)  +  AU  (r)  (III. 14) 

D  o 

where  U  (r)  is  the  usual  spherical  optical  potential  and  AU  is  the  non- 
o 

spherical  part  of  U^.  It  is  clear  from  equation  (III. 12)  that  AU  now 
becomes  the  transition  operator  for  the  excitation  of  the  state  We 

can  rewrite  the  T-matrix  element  (equation  (III. 12))  for  a  single  exci¬ 
tation  of  2^-pole  collective  surface  mode  of  an  even  target  nucleus  as 
(Au65  and  Sa65) : 

Tfi  =  <X(-} (kf ,r) | <IM| AU | 00> | X (+) (k.,r) >  (III. 15) 

where  |00>  describes  the  ground  state  with  spin  parity  0+.  | IM>  des¬ 

cribes  a  collective  excited  state  with  spin  I  and  projection  M.  The 
parity  of  the  state  is  given  by 

1Tf  =  (-)1.  (III. 16) 


In  this  chapter  there  is  necessarily  a  lot  of  algebra.  Not  all  the 
steps  are  included,  however,  the  interested  reader  may  find  some  of  the 
missing  details  in  the  thesis  of  Sherif  (Sh68) .  It  must  be  noted  that 
because  of  the  spin-orbit  interaction  in  equation  (II. 1),  the  x's  are 
matirces  in  the  particle  spin  space.  For  example 


X(+)(k.,?)  =  l  I \  y  '>X^  (!>?)  (III. 17) 

u[  1 

where  y^  is  the  spin  projection  of  the  incident  particle  and  |y  y^>  is 
the  spin  eigenfunction.  The  quantization  axis  is  taken  to  be  along 
the  incident  beam  direction.  Terms  with  y^  ^  }_u  allow  the  possibility 
of  spin  flip  during  the  elastic  scattering.  The  same  thing  holds  for 
^ .  One  can  write  the  matrix  element  as 

Tfi=,,£  <Xy^yf(^f>A  |<J  uy<IM|AU  I  00>|y  u^>|X^  (!,?)>  (HI. 18) 


UfUf 
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where  is  the  spin  projection  of  the  scattered  particle.  The  XTs  are 
distorted  waves  in  the  relative  coordinate  of  the  projectile  and  the 
target.  They  are  solutions  of  the  Schrodinger  equation,  with  appropri¬ 
ate  boundary  conditions,  for  a  spherical  optical  model  potential. 

The  interaction  potential  AU  is  capable  of  connecting  the  ground 
state  to  the  collective  states,  and  usually  is  obtained  by  deforming 
the  optical  potential  UQ(r).  This  could  be  achieved  by  taking  the 
density  function  fCtjR^a^)  (see  equation  (11.28))  and  allowing  its 
radial  parametrs  to  depend  on  the  nuclear  orientation  such  that 

Rf  +  R  +  a(r)  (III. 19) 


where  R  is  the  radius  parameter  of  a  typical  spherical  optical  poten¬ 
tial,  and  a(r)  is  the  displacement  of  the  nuclear  surface  in  the  direc¬ 
tion  r  and  may  be  given  the  multipole  expansion  (Au65) : 


*M 


«(»-  l 

L,M  L,N  L 


(III. 20) 


where  £T  is  the  dynamical  coordinate  which  describes  the  surface  dis- 
L,M 

placement.  If  we  expand  f  (r  ,R +a(r“)  ,a)  and  keep  terms  to  first  order 
in  a(£)  (B164) ,  we  get 


f  (r  ,R +a(r)  ,a)  »  f(r,R,a)  +  a(r)  ||-  . 


It  is  convenient  to  write 


(III. 21) 


AU  =  AU_  +  AU  +  AU.  +  AU  (III. 22) 
C  r  i  SO 

referring  to  Coulomb,  real,  imaginary  and  spin  dependent  parts,  res¬ 
pectively.  The  contribution  to  AU  in  terms  of  the  optical  potential 
parts  are  (see  equation  (11.27)) 


AUr  =  -VQao(e)  gf-  f(r,Ro,ao) 
o 


(III. 23) 
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AUi  =  ^.(r)(W-4a.WD  JL  f(r,R.,a.) 

1 


(III .24) 


AU 


SO  m 


1TC- 


tvso°-  v 


-  ... 
aso(r)  - 


3R 


so 


x4  V 
1 


+  iwsoa-  7l“soi(r) 


3£(r’RSOI>V 


3R 


SOI 


1  -*■ 
V] 

1 


(III. 25) 


where  AU  comes  from  the  deformation  of  the  spin  dependent  potential, 
SO 

it  is  of  the  full  Thomas  form  (Sh69) ,  and  it  could  be  split  into  two 
parts  as  follows 


AUS0  ^s,def(l)  +^s,def(2) 


(III. 26) 


where 


U 


s ,def (1) 


ft 


TTC' 


2  rTT  1  9  I  3f 

[VS0aS0(r)  r  3r  [sr^ 
3f 


lWSOiaSOI(r)  r  3r 


3r 


SOL 


a  •  L 


(III. 27) 


U 


s,def (2) 


m 


of 


TT 

+  iW 


2  [vso  l!5-J-^“so(e»xT^ 

bU 


(III. 28) 


It  is  found  that  the  contribution  to  the  deformed  potential  due  to  AU^ 
is  important  in  this  calculation.  It  improves  the  calculated  cross- 
sections  at  backward  angles,  and  also  it  is  necessary  to  produce  the 
inelastic  polarization  (Sh69) .  The  Coulomb  contribution  to  AU  is  given 

by 


<im[auc|oo> 


3ZZ’e2 


r(r/R  )  r  <  R 

sm  Y*V)I  c 

(21+1)  I()l(R/r)I+l  r>R, 


(III.  29) 


18 


where  Z  and  Zf  refer  to  the  atomic  numbers  of  projectile  and  target 
respectively,  and  the  spherical  charge  distribution  is  deformed  by  let¬ 
ting  R  -*R  +  a  (r),  but  maintaining  the  uniformity  of  the  distribution, 
c  c  c 

The  contribution  from  the  inner  region  (r  < R^)  to  equation  (III. 29)  is 
found  to  be  very  small. 

The  nuclear  structure  information  is  contained  in  the  matrix  element 


«&  Jau|$.>. 

f 1  1  i 


(III. 30) 


If  we  consider  the  exictation  of  a  rotational  level  in  an  even-even 
nucleus,  i.e.  ( = 0+) ,  then  the  wave  functions  are  given  by 


l 


00>  = 


IM> 


Btt 
21+1 


1/2  hntD>> 


Btt" 


1/2  *■ 

rint  0M 


(111. 31) 

(111. 32) 


where  D  is  a  rotation  matrix,  e  denotes  the  orientation  of  the  body 

axes  w.r.t.  the  space  fixed  axes  and  ib .  is  the  intrinsic  wavefunction. 

hnt 

Invoking  axial  symmetry  one  can  rewrite  equation  (III. 20)  in  the  fol¬ 
lowing  fashion 


(III. 33) 


where  is  real  and  normally  written  as 
deformation  parameter  and  (0f,cj)f)  refer  to  body  fixed  system.  Then  in 
the  space  fixed  system  one  can  write 


Y°(e’,r)  =  l  4i<e)Y*M0f) 


(III. 34) 


and 
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1/2 

^flaJV  ■  (2^2  ^intDOMf(£)lX  Ro6LDOM(£) 
8tt  f  LM 


I  Ro6jfrr) 


8tt  LM 


<D0Mf(e>lD0M(£)lD00(£)>  • 


Since,  (Ro57) 

<D0M  <e)|  DOM(£)l  D00(e)>  =  ^r6L,I6M;1 


(III. 35) 


(III. 36) 


then 


R  3 

<®  let  |®.> - °  -  ■  Y*M(r)6  6 

f  o'  1  (2I+1)1'2  L  I,L  Mf,M 


*M 

c(i)YiM(e) 


(III. 37) 


Now,  if  we  define  a  radial  form  factor: 
R 


F  (r)  =  f  (r,R  )  -  V  f(r,R  ,a  ) 

o  R  c  c  o  3R  o’  o 

o  o 


-  i(VV(w-4aiwD^ 

i 


(III. 38) 


where  fc(r,R^)  is  obtained  from  equation  (III. 29)  by  omitting  the  part 
*M 

C(I)Y  ,  then  for  the  Coulomb,  real  and  imaginary  contributions  one  gets 


I  ,  *M 

<IM|AUc  +  AUr  +  AUjOO  =  C(I)Fo(r)YI  (r) 


(III. 39) 


Similarly,  for  the  spin-dependent  parts: 


<D1lUs,def(l)l°0>  "  CS0(I)Fl(r)YlM^)?-J 


(III. 40) 
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:IM|U 


s ,def (2) 


1  -*■ 

|00>  =  Cg0(I)F2(r)  O'  [V(YI  )  x-V] 


(III. 41) 


Fl(r)  Im  cj  [VSO  r  3r3R  f (r,ESO’aSO) 
7T  J  SO 


+  iwsoi  r  57al^7f(r’Rsoi>asoI)]  (III‘42) 


F2(r)  = 


m  c 
TT 


Xo  wr  f(r’Rso’aso)  +iWsoi  f(r>Rsorasoi)1-(III-43) 


SO 


SOI 


The  T-matrix  element  may  now  be  written  as: 


fi 


=  + 


Tl  +  T2 


where 


T0  =  C(I)  l  ,  1 ^  u^|Fo(r)Y*Mce)|i  ui>; 


f  f 


*l.vki,r)> 

ii 


(III. 44) 


Ti  =  cso(I)  l,  <*2  <£ f,r)|<| 


u:Uf 


^11 


(III. 45) 


T., 


CSO(I)  l. 


<x.(T) 


UfUf(Sf,?)i4u-|F2wa-[v(Y*M)4v] 


(III. 46) 


These  expressions  are  also  valid  for  the  excitation  of  a  vibrational 
state  | IM>. 

->■ 

If  we  choose  the  z-axis  along  the  vector  (see  Figure  1) ,  then 
one  can  write  the  following  partial  wave  expansion  of  the  distorted 


. 
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Figure  1.  Coordinate  system. 


waves 


yiyi  1  ki  ij 


l  i^e  Sl(2l+l)1/2  <Jl|-;Oyi|jui> 


..1  ...  FJM(ki’r) 


(III. 47) 


and 


*(_\  ->  +  _  _o  »  ia(?  » 

xu’u>(kf,r)  =r  J  1  6  V  f(0’o) 

ufuf  *  kf  Jl'j'v  Z 


<i'  f,v-yf  y f | J ’ vxfi, '  y;v-u'  yj|j'v> 


*v-y*  ^  F.,., (k  ,r) 

X  Y  ,  f(r)  *  3  r  f 


(III. 48) 


Where  9  is  the  scattering  angle,  <&s,m  m  |jm.>  denotes  the  Clebsh- 

36  s  J 

Gordan  coefficients,  a.  is  the  Coulomb  phase  shift 


On  =  arg  T(Z  +  1  +in) 


(III. 49) 


where 


ZZ'e2 

h2k 


(III. 49) 


and  the  F^_.  are  regular  radial  wave  functions  and  are  normalized  such  that 
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^F*j=i[H  -V*1 


(III. 51) 


where  ri0  .  are  the  reflection  coefficients  or  scattering  matrix  elements 
for  the  (ilj)  partial  wave,  andH^  is  the  outgoing  wave  Coulomb  function  (Hu59) 
The  final  expression  for  the  T-matrix  is  given  by  equation  (36) 
in  reference  (Sh69) .  In  terms  of  the  T-matrix  element,  which  is  writ¬ 


ten  as  T(yi,]if  ,M) ,  we  can  write  the  differential  cross-section  in  the 


following  fashion 
*2  k 


da 

dft 


f  1 


2  |TCu±;ufM)| 

4tt  fi  i  y.y  M 
if 


(III. 52) 


where  m  is  the  reduced  mass,  and  the  factor  —  arises  from  averaging 
over  initial  spins.  The  summation  over  y^  in  equation  (III. 52)  can  be 


taking  care  of  using  the  relation  (Sh71) 
l  |T(ui;UfM)  I2  =2  l  |  T(y,UfM)  | 


(III. 53) 


UiSUfM 


UfM 


The  inelastic  polarization  could  be  written  also  in  terms  of  the 
T-matrix  element.  But  before  doing  that  let  us  define  a  density  matrix 
p  in  the  following  fashion  (W056) 


p  =  TT 


(III. 55) 


which  is  a  matrix  in  the  spin  space  of  the  scattered  particle.  If  we 


write  T(y.;y  M)  as  a  matrix  element  T  (M)  such  that 
1  t  yfUi 

(T(M)T+(M))UfU,  -  l  T  (M)T*  00 

Mi 


(III. 56) 


and  making  an  incoherent  sum  over  M,  we  can  write  the  following 


P,.  ..»  =  l  T(y  ;u  M)T  (y  ;y'M)  . 
MfMf  y.M 


(III. 57) 


The  polarization  vector  P(0)  is  the  expectation  value  of  the  a  operator 


P(0) 


trpg  _  tr  TT  a 
trp 


(III. 58) 


tr  TT 

where  the  denominator  is  given  by 


trp  =  tr  TT 


=  l  |t(u. ;ufM) | 2  (III. 59) 

U±VfM 

and  is  thus  proportional  to  the  differential  cross-section,  and  the 
final  expression  for  the  polarization  could  be  written  as 

P (9)  =  - 2 — .  (III. 60) 

I  |x(|;wfH)|2 

UfM 


It  was  found  that  this  model  is  successful  in  explaining  the  experimen¬ 
tal  data  for  low-lying  excited  collective  states  at  low  energies.  (Ba- 
62,  Fr65 ,  Fr67,  Fu68,  Sh69,  Wi68) .  In  the  present  work  we  have  calcu¬ 
lated  both  the  inelastic  cross-section  and  polarization  for  different 
targets  at  intermediate  energies.  The  results  of  these  calculations 
according  to  this  model  will  be  shown  and  compared  with  the  results  of 
the  unconventional  potential  calculations  in  Chapter  VII. 


. 


IV.  DIRAC  PHENOMENOLOGY  FOR  NUCLEON  ELASTIC  SCATTERING 


We  have  seen  in  Chapter  II  how  to  do  elastic  scattering  according 
to  a  non-relativistic  model,  i.e.  the  standard  optical  model.  In  this 
chapter  we  discuss  another  model  for  doing  elastic  scattering  also,  but 
from  a  different  point  of  view.  A  relativistic  optical  model  is  con¬ 
structed  to  be  used  for  the  analysis  of  proton  elastic  scattering  (Ar81) . 
The  model  is  based  on  the  Dirac  equation  with  a  mixture  of  Lorentz  sca¬ 
lar  potential,  Ug,  and  the  time  like  component  of  a  four-vector  poten¬ 
tial,  U^.  The  time  independent  Dirac  equation  which  describes  proton 
elastic  scattering  is  (ft  -  c  ~ 1) 


{a-p  +  B(m+Us(r))  +  (Uv(r)  +Vc(r)  h|j(r)  =  Et()(r) 


(IV. 1) 


where  a  and  3  are  the  4  by  4  Dirac  matrices  and 

3  = 


0  6' 
a  o  ’ 


ri  o 

o  -i 


(iv. 2) 


a  is  the  Pauli  2  by  2  spin  matrix,  I  is  the  2  by  2  unit  matrix,  p  is  the 
momentum  operator  of  the  incident  proton,  m  is  the  rest  mass  of  the 
proton,  i(j  is  the  Dirac  spinor,  V^  is  the  Coulomb  potential  (it  is  de¬ 
termined  from  the  empirical  nuclear  charge  distribution  given  by  elec¬ 
tron  scattering  experiments)  and  E  is  the  proton  total  energy. 

Major  properties  of  the  nucleon-nucleon  interaction  can  be  repre¬ 
sented  by  one  boson  exchange  potentials  employing  the  known  mesons  (tt, 

TT  *4* 

p  ,co,cj) , . . . )  plus  the  addition  of  the  J  =0  ,  T=0  "a"  meson  which  is  presently 
interpreted  in  terms  of  two-pion  exchange  processes.  The  dominant  charac¬ 
teristics  of  the  nucleon-nucleus  interaction  for  a  spin  zero  isospin 
zero  target  nucleus  are  expected  to  be  represented  by  the  exchange  of 
neutral  scalar  and  vector  mesons  (Ar81) . 


24 


' 


' 


25 


The  main  effect  of  using  a  mixture  of  a  Lorentz  scalar  and  the  time  like 
component  of  a  four-vector  potential  in  the  Dirac  equation  was  original¬ 
ly  discussed  by  Furry  (Fu36) .  He  pointed  out  that  to  lowest  order  in 
v/c  the  sum  uv  +ug  contributes  to  the  central  potential,  while  the  dif¬ 
ference  -Ug  appears  in  the  spin  orbit  potential.  This  feature  of 
the  potential  mixtures  plays  an  essential  role  in  the  description  of  both 
elastic  scattering  and  polarization  data  at  intermediate  energies  (Ar81) . 
This  model  is  quite  successful  in  describing  the  experimental  data  at 
intermediate  energy.  In  particular,  for  energies  above  200  MeV,  it 
seems  to  be  superior  than  the  standard  optical  model. 

In  order  to  compare  with  the  non-rel-ativistic  optical  model,  it  is 
convenient  to  rewrite  equation  (IV. 1)  in  second  order  form.  If  we 
write  the  Dirac  spinor  in  terms  of  an  upper  and  a  lower  component  as 
follows 


and  if  we  substitute  equations  (IV. 2,  IV. 3)  into  equation  (IV. 1),  we 
get- 


7o  ct-p 

fm+U 

s 

0 

ru  +v 

v  c 

0 

V 

'  v 

.la*P  0 

+ 

o 

-(m+Us). 

+ 

0 

u  +v 

V  C‘ 

w 

=E 

.  H 

(IV. 4) 

This  will  lead  to  two  coupled  equations  as  follows 

a-p^  +  (m+Us+Uv+Vc)^  =  (IV. 5) 

+  «JV  +Vc  -m  -Us)!^  =  . 


(IV. 6) 
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From  equation  (IV. 6)  we  get 


£ 


_ 1 _ 

E+m+U  -U  -V 

SVC 


a-p  iD 


1 

D 


a-p  \p 


y 


(IV.  7) 


where 


DEE+m  +  U  -  U  -V  . 

SVC 


(IV.  8) 


Using  equation  (IV. 7)  in  equation  (IV. 5)  we  get 


+  (m+Us+Uv+V’i'u  =  *V 


(IV. 9) 


After  some  mathematical  manipulations  one  can  write  the  following 


a-p  -  a.pij^  =  (a-p  -)  a-pip  +  -  (a-p)  (a-p)ijj 

i  1  3D  ■*  1  1  3D 

=  77  3?r-P,J'u  -77  37  a-L’*'y 

+  i  P2*  . 

D  H  ry 


(IV. 10) 


Substituting  equation  (IV. 10)  in  equation  (IV. 9),  one  can  rewrite 
the  latter  as  follows 

P2^  +D(m+U  +  U  +V  )i|j  -  ^--^-a-L^.. 

Tu  svc  Ty  D  r  dr  y 


D  r  3r  ry 


(IV.  11) 


To  eliminate  the  Darwin  term  —■  —  ■  r»pii  one  can  make  the  substitution 

D  r  3r  ry 


1  /  0 

t(r)  =  (D(r) )  $  (J)  . 

y  V 


(IV. 12) 


Then  the  term  P  ip  could  be  written  as 


p2ip  =  p2d1/2$ 

y  y 


(p 


2  1/2,,  ,  ->1/2  +  ,  1/2  2 

D  )  $  +  2PD  *P<1>  +  D  P  $ 


(IV. 13) 
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and  after  some  algebra  it  could  be  rewritten  as 

p2*  =  -iD-1/2D'*  +-U-3/V2*  -^d-1/2d"$ 

y  r  y  4  y  2  y 


_d-1/2d.^+d1/2  2$ 

9r  y 


(IV. 14) 


where  (f)  and  (")  mean  the  first  and  the  second  derivatives  with  respect 
to  r,  respectively.  Equation  (IV. 11)  could  be  rewritten  as 

P2*  (r)  +  J  (D’/D)2$u  -  t(D"/D)$u(r)  -  t(D •  /D)<yr) 

+  D(m+Us +Uv+Vc)$iJ(r)  - t(D' /D)a-L$u(r)  =  ED^ (r) .  (IV .15) 

Since 


D(m+U  +U  +V  )  -  ED  =  2E(U+V  )  +  m2  -  E2 +2mU 
SVC  c  s 

+  u2  -  (u  +v  )2  . 

S  VC 


(IV. 16) 


Thai  equation  (IV. 15)  reduces  to 

(P2  +  t*  (D'/D)2  (D”/D)  --  (D'/D)  +2mU  +  2E(U  +V  ) 

4  2  r  svc 

+  U2  -  (U  +V  )2  --  (D'/D)a.L>$  (r)  =  (E2  -m2)$  (r)  .  (IV. 17) 
s  v  c  r  y  y 

If  we  divide  equation  (IV. 17)  by  2m,  noting  that 


2  2 
B  -m 

2m 


k__ 

2m 


(IV.  18) 


one  can  rewrite  (IV .17)  as  follows 

2  .  . 


(k  +Ua,yr>  +U5n(r)3-L)0,1(r) 


ef  f 


SO 


u 


k/ 

2m  ‘y 


(IV. 19) 


where  ue£f(r)  is  the  effective  central  potential  defined  as  follows 

+  ^-  (u2  -  (u  +v  )2] 

2m  s  v  c 


and  Ug^(r)  is  the  spin  orbit  potential  defined  as  follows 


(IV. 20) 


‘■'vA' 
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tt  r  \  M  1  dD 

US0<-r)  2m  D  r  dr  • 


(IV. 21) 


Equation  (IV .19)  looks  like  the  non-relativistic  Schrodinger  equation 
which  is  usually  written  in  the  following  fashion 
2 

if;  +  Vijj  =  El{J  .  (VI. 22) 

zm 


Since  equations  (IV. 19)  and  (IV. 22)  produce  equivalent  elastic  scat¬ 
tering,  equation  (IV. 19)  is  called  the  Schrodinger  equivalent  equation. 
It  should  be  mentioned  here  that  the  wave  function  ,  even  though  it  is 
not  the  upper  component  of  the  Dirac  spinor,  it  has  the  same  asymptotic 
behaviour  as  i|j  in  the  absence  of  the  Coulomb  potential,  this  is  because 

of  the  short  range  of  U  and  U  . 

vs 

Equations  (IV. 20,  IV. 21)  show  that  if  is  repulsive  (as  it  would 
be  if  it  is  due  to  neutral  vector  meson  exchange) ,  and  is  attractive 
(as  it  would  be  if  its  origin  is  neutral  scalar  meson  exchange) ,  then 
and  Ug  tend  to  cancel  in  the  effective  central  potential  while  they 
add  in  the  spin-orbit.  Notice  that  there  is  no  explicit  spin  depen¬ 
dence  in  the  optical  potentials  and  U  ,  where  those  are  usually  writ¬ 
ten  as 


U  (r)  =  V  f(r,R  ,a  )  +  iW  f (r,R  ,a  ) 
v  v  v  v  v  wv  wv 


(IV. 23) 


and 


U  (r)  =  V  f(r,R  ,a  )  +  iW  f(r,R  ,a  ) 
s  s  s’  s  s  ’  ws’  ws 


(IV. 24) 


where  the  function  f  (r ,R  , a  )  is  of  the  Woods-Saxon  form.  Note  that  the 

x  x 


effective  spin-orbit  potential  is  completely  specified  by  the  potential  mix¬ 
ture.  We  will  call  this  type  of  calculation  "the  Dirac  model"  in  the 
present  work.  The  characteristics  of  the  potentials  and  the  elastic 


■ 


' 
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scattering  cross-section  and  polarization  predicted  by  this  model  will 
be  discussed  and  shown  in  Chapter  VI,  and  they  will  be  compared  with 
their  counterparts  from  the  standard  model. 


V.  RELATIVISTIC  INELASTIC  SCATTERING  (THE  DIRAC  MODEL) 


We  have  seen  in  chapter  III  how  to  do  inleastic  scattering  accord¬ 
ing  to  a  non-relativistic  model,  namely  the  extended  standard  optical 
model.  This  model  has  been  successful  in  describing  the  experimental 
data  at  low  energies,  but  its  credibility  becomes  suspect  at  intermedi¬ 
ate  and  high  energies.  The  intermediate  energy  range  attracts  the  in¬ 
terest  of  many  nuclear  physicists  these  days,  especially  with  the  avail¬ 
ability  of  facilities  like  SACLAY,  SIN,  LAME ,  IUCF,  and  TRIUMF .  Since 
the  Dirac  equation  has  been  successful  in  describing  proton  elastic 
scattering  at  intermediate  energies  in  recent  years  (Ar81,82),  then  the 
question  arises  as  to  whether  there  are  serious  implications  for  in¬ 
elastic  scattering  as  a  result  of  the  unconventional  shape  of  the  real 
central  potential.  Satchler  (Sa83)  has  recently  studied  the  effect  on 
the  inelastic  cross-section  for  some  low-lying  states  near  200  MeV,  but 
his  calculations  were  based  on  varying  the  shapes  of  the  potentials  and 
using  the  non-relativistic  Schrodinger  equation.  In  this  chapter  we 
will  discuss  a  model  based  on  the  use  of  the  Dirac  equation  to  describe 
elastic  scattering  as  discussed  in  chapter  IV.  To  allow  for  inelastic 
scattering,  we  have  to  deform  the  scalar  and  the  vector  potentials 
for  this  purpose  we  use  the  same  method  used  earlier  to  deform  the 
non-relativistic  potentials  (see  chapter  III) .  The  vector  and  scalar 
potentils  are  defined  by  equations  (IV. 23)  and  (IV. 24)  respectively. 

We  can  write  the  deformed  vector  potential  in  the  following  fashion 

UD(r)  =  V  f(r,R  +a  (r)  ,a  )  +  iW  f(r,R  +a  (r),a  )  (V.l) 

v  v  v  v  v  v  wvwv  wv 

and  again  express  this  in  terms  of  spherical  and  non-spherical  parts: 

n®(r)  =  Uv(r)  +  AUv(?)  (V.2) 
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if  we  expand  the  function  f  and  keep  terms  of  first  order  in  ot(r)  we 
can  write  (see  equation  (III. 21)) 


AU  (r)  =  a  (r)V  +  iW  (r)W 

v  v  v  R  v  v  3R 

v  wv 


(V .  3) 


and  by  the  same  method,  we  can  write  for  the  scalar  potential 

AU  (r)  =  a  (r)V  -  +  ia  (r)W  . 

s  s  s  3R  ws  s  3r 

s  ws 


(V  .4) 


The  effective  central  potential  is  defined  by  equation  (IV. 20),  and  it 
could  be  split  into  two  parts  as  follows 


Ueff(r)  =  A(r)  +  B(r) 


(V  .5) 


where 


,  (he)  r3  ,D',2  1  ,D",  1  ,D' 

A(r)  =  “ST  [4  (T}  -  2  -  r  (T>] 


(V  .6) 


and 


B(r)  =  U  (r)  +  |  (U  +V  )  +-±-  [U2  -  (U  +V  )2] 
s  mvc2ms  vc 


(V.7) 


where  D  is  defined  by  equation  (IV. 8).  The  deformation  of  B(r)  follows 
as  a  result  of  deforming  the  scalar  and  vector  potential  and  it  could 
be  written  as  follows 

BDdo  =  U°(r)  +  -  (UD(?)  +  V°(r))  +  J-  [U  °(.r) 
s  m  v  c  2m  s 

2n  D  ->•  D  2r\  -> 

-  U  U(r)  -  2UU(r)VU(r)  -  V  U(r)]  .  (V.8 

V  VC  c 

which  may  be  written  as 

BD(r)  =  B(r)  +  AB(r)  (V.9) 


where 
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AB(r)  =  AU  (r)  +  -  [AU  (?)  +AV  (?)  J 
s  m  v  c 

+  [2U  (r) AU  (?)  -  2U  (r)AU  (?)  -  2U  (r) 

2m  s  s  v  v  v 

AV  (?)  -  2V  (r)AV  (?)  ]  .  (V.10) 

c  c  c 

The  contributions  to  AB(r)  due  to  Coulomb  effects  were  investigated  and 
found  to  be  negligible,  so  we  decided  to  leave  these  terms  out,  this 
will  be  discussed  in  chapter  VII.  We  can  then  rewrite  equation  (V.10) 
as  follows 

AB(r)  =  AU  (?)  +  -  AU  (?)  +  -  [U  (r)AU  (?) 

s  m  v  ms  s 


-  Uv(r)AUv(r)] 


(V.ll) 


The  deformation  of  A(r),  depends  on  the  deformation  of  D(r)  and  its 
derivatives,  and  it  could  be  done  in  the  same  way  mentioned  before. 

One  can  write 


AD(r)  =  AU  (r)  -  AU  (r) 
s  v 


AD'(r)  =  AUT(r)  -  AU  1  (r) 
s  v 


(V .  12  ) 

(V .  13) 


and 


AD"(r)  -  AU"  (r)  -  AU"(r) 
s  v 


(V .14) 


After  some  straightforward  algebra  one  can  write  AA(r)  as  follows 


4A(?)  -  nib  [f  (t)AD'  -  f  + 1 

-h*f>  -?(T>  +7d,(7>] 

if  we  use  the  following  relation 

V2D  =  -  Df  +  D" 
r  r 


(V  .15) 


(V .  16) 
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then  we  can  write 


and 


D"  =  V2D  -  -  D1 
r  r 


(AD)"  =  V^(AD)  -  |  (AD)' 


(V  .17) 


(V  .18) 


Substituting  equations  (V.17),  (V.18)  in  equation  (V.15)  yields 


AAW  (^-Had)'  -f  02(^)] 


+  0)£)_  [V^D(^j)  -  £  V^(AD)] 


4m 


(V .19) 


Then  the  deformed  piece  of  the  effective  central  potential  AU^^Cr) 
just  the  sum  of  equations  (V.ll)  and  (V.19)  since 


Au  .-(r)  =  AA(r)  +  AB(r) 
eft 


(V  .20) 


If  we  define  a  formfactor 

3f  3f  i  ?  3f  o  3f 

Vr)=Vs)T+Vvjf+i  [VS  fs  Sir  -  Vws  3R^ 

S  V  s  ws 


3f  -i  3f 

+  i[W  +  »  (f  — — 

s  3R  m  s  s  s  3R 


3f 


ws 


+  fws  ar)] 
ws  s 


(he) 


D\  3 
D 


D’D"  ,  1  „2 
D 


+  —  [3 (— )  -  3  ^-5-  +  V"”  (D 1 )  -  V"D] 


D_|_  „2, 
2 


(V .21) 


then 


<IM|AU  (r) | 00>  =  C(I)F  (r)Y_  (r)  . 


(V .22) 


eff 1  ”  ~ o  ~I 

The  deformation  of  the  spin  dependent  part  is  more  involved  than 
that  of  the  effective  central  potential  but  it  is  done  by  the  same 
method.  The  deformed  spin-orbit  potential  is  given  by 


’“so®  ■-  -V 5-w»  *; 

D(r) 


(V.23) 
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where 


D(r) 


AP(r) 

D2(r) 


(V.24) 


D(r)  D(r)+AD(r) 
if  we  substitute  equation  (V.24)  into  equation  (V .23)  ,  this  yields 
2 


<,*>  ■  ^  >< 1 


^rr-r  -  V(D(r)  +  AD(?))  xp] 

UW  D  (r) 


(V  .25) 


which  could  be  rewritten  as  follows 


-  a-7(D+AD)  xp  ] 

D 


(V  .26) 


The  first  term  on  the  r.h.s,  is  the  spherical  spin-orbit  interaction,  and 
the  other  parts  form  the  deformed  spiece. 
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where  the  last  term  drops  out  if  we  keep  AU^Cr)  to  first  order  in  a(£)» 
and  the  second  term  is  not  hermetian  when  the  potentials  are  real  and 

we  have  found  its  effect  on  the  calculations  to  be  negligible  since  it 

2 

is  divided  by  D  ,  so  we  decided  to  leave  this  term  out.  Then  equation 
(V.27)  could  be  rewritten  as  follows 
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where 


U 


s , def (1)  2m 


o  2 

i  3  f  3  f 

-•(he;  1  r/  s  v 

— t*-  tr  l  (a  (r)V  r--r-_  -  a  (r)V  r— — 

D  s  s  3rdR  v  v  9r3R 


+  i(aws(?)Ws  Sf- 
ws 


32f 


-  a  (r)W 

wv  v  3r3R 


-)  ]a-L 


(V  .29) 


and 
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which  is  of  the  full  Thomas  form  (Sh68) .  We  can  define  the  form- 


factors  for  U  ,  and  U  ,  as  follows 
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Then  if  we  substitute  equations  (V.21),  (V.33)  and  (V.34)  in  equations 
(III. 44),  (III. 45)  and  (III. 46)  the  calculations  will  follow  the  same 
order  as  that  of  chapter  III. 

We  have  done  the  calculations  for  the  inelastic  cross-section  and 
polarization  according  to  this  model.  We  will  compare  the  results  of 
this  model  with  those  obtained  from  the  non-relativistic  calculations 


which  are  based  on  the  extended  standard  optical  model,  and  with  the 
results  of  another  model  (Sa83)  for  some  cases  in  chapter  VII. 


VI.  ELASTIC  SCATTERING  FITS 


In  this  chapter  we  shall  discuss  the  elastic  scattering  fits  for 
different  even-even  nuclei  at  different  bombarding  energies.  Two  dif¬ 
ferent  codes  were  used  in  fitting  the  elastic  cross-section  and  polari¬ 
zation  data.  The  first  code  is  a  standard  optical  model  potential 
called  Magali  (Ra69) ,  relativistic  kinematics  were  used  in  this  code, 
i.e.,  we  used  an  effective  bombarding  energy  that  ensured  the  relati- 
vistically  correct  centre  of  mass  momentum.  The  Dirac  equation  based 
potential  has  been  taken  from  the  analysis  of  Cooper  et  al.  (Co82) . 

This  analysis  was  done  using  the  code  Runt  (1981) .  The  code  Magali 
is  based  on  the  optical  potential  described  by  equation  (11.27)  while 
the  code  Runt  is  based  on  the  optical  potential  described  by  equations 
(IV. 20)  and  IV. 21).  In  both  codes  12  parameters  were  used  in  an  auto¬ 
matic  search  routine  using  a  chi-squared  minimization  technique 
(see  equations  (II. 30)  and  (11.31)).  These  parameters  were  varied  to 
fit  the  calculated  proton  elastic  scattering  cross-section  and  polari¬ 
zation  to  the  experimental  data.  The  radius  of  the  charge  distribution 
R^  was  held  constant  in  both  codes,  since  the  fit  is  fairly  insensitive 
to  it.  The  parameters  obtained  from  the  code  Magali  are  listed  in 
Table  I  for  different  cases.  We  show  a  sample  of  the  parameters  ob¬ 
tained  from  the  code  Runt  (Co82)  in  Table  II.  In  both  tables  we  have 
listed  the  optical  potential  parameters  for  different  target  nuclei  at 
different  energies.  As  can  be  seen  from  Table  I  the  parameters  of  the 
potential  are  energy  dependent,  they  vary  with  increasing  energy.  They 
also  depend  on  the  target  mass  number  as  it  was  reported  by  Nadaseen 
et  al.  (Na81)  and  by  van  Oers  et  al.  (Va74) .  We  also  show  the  calcu¬ 
lated  reaction  cross-section  (cO  in  millibarns  for  each  case. 
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TABLE  I 


OPTICAL  MODEL  PARAMETERS 
(THE  STANDARD  MODEL) 


Target 

12c 

24Mg 

28„. 

Si 

32s 

40„ 

Cel 

4CL 

Ca 

40^ 

Ca 

40 

Ca 

208Pb 

E  (MeV) 
P 

200 

155 

155 

155 

155 

181 

300 

500 

200 

V  (MeV) 
o 

14.18 

12.29 

13.96 

12.90 

13.68 

17.52 

25.72 

-41.69 

13.17 

ro(F) 

1.30 

1.39 

1.34 

1.37 

1.37 

1.29 

1.07 

0.70 

1.33 

ao(F) 

0.68 

0.57 

0.50 

0.59 

0.61 

0.70 

0.68 

0.79 

0.58 

W(MeV) 

17.64 

22.31 

58.56 

19.55 

6.00 

13.81 

20.26 

54.21 

25.90 

WD(MeV) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

rw(F) 

0.95 

0.98 

0.78 

1.06 

1.26 

1.15 

1.19 

1.04 

1.13 

a  (F) 
w 

0.78 

0.81 

0.78 

0.80 

0.60 

0.81 

0.78 

0.79 

0.88 

VS0(MeV) 

4.01 

4.19 

14.1 

4.72 

4.86 

3.20 

3.58 

4.95 

3.38 

vso(F) 

0.88 

0.90 

0.42 

0.93 

1.06 

1.02 

1.06 

0.99 

1.10 

aso(F) 

0.66 

0.71 

0.96 

0.70 

0.58 

0.65 

0.62 

0.75 

0.73 

wsoi 

-2.09 

-1.65 

-5.94 

-1.24 

-1.61 

-2.25 

-1.39 

-4.15 

-2.69 

(MeV) 

rsoi(F) 

0.96 

0.96 

0.78 

1.02 

1.06 

1.01 

0.94 

1.08 

1.08 

tn 

O 

H 

/■ — s 

0.37 

0.57 

0.59 

0.58 

0.68 

0.56 

0.60 

0.68 

0.70 

V® 

194.74 

402.5 

501.2 

494.8 

296.8 

495.5 

568.5 

591.5 

1827.7 
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TABLE  II 
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OPTICAL  MODEL  PARAMETERS  FOR  Ca  AT  E  =181  MeV 

P 

(THE  DIRAC  MODEL) 


— 

V^(MeV) 

ry(F) 

av(F) 

295.184 

1.0416 

0.6017 

W  (MeV) 

r  (F) 
wv 

a  (F) 
wv 

-89.7876 

1.0774 

0.6518 

V  (MeV) 
s 

rs(F) 

as(F) 

-404.386 

1.0288 

0.6314 

W  (MeV) 
s 

r  (F) 
ws 

a  (F) 

WS 

93.6447 

1.0871 

0.6262 

a  (mB)  =  594.80 
R 
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In  the  present  work  we  will  label  the  curves  resulting  from  the  stan¬ 
dard  model  as,  STANDARD,  while  those  resulting  from  the  Dirac  model  as, 
DIRAC.  The  curve  labeled  Dirac  is  actually  obtained  using  the  Schrodinger 
equivalent  potential  (the  effective  central  potential  (r)  ,  see  equa¬ 
tion  (IV. 20))  as  well  as  relativistic  kinematics  in  the  Schrodinger  equa¬ 
tion.  The  wave  function  which  satisfies  the  Schrodinger  equation,  with 
Ueff  (r)  ,  differs  from  the  upper  component  of  the  Dirac  wave  function  only 
in  the  nuclear  interior  but  they  have  the  same  asymptotic  behaviour .  One, 
therefore,  expects  the  equivalent  potential  (U^^  (r) )  when  used  in  a 
Schrodinger  equation  based  elastic  scattering  calculation,  to  reproduce 
the  same  elastic  scattering  obtained  using  the  Dirac  equation.  This 

point  has  been  checked  by  actual  numerical  calculations  for  targets  up 
208 

to  Pb  by  Cooper  (Co81) . 

We  will  show  below  a  few  cases  from  the  analysis  of  proton  elastic 
scattering  in  this  work. 

vi. i  12c 

Figures  2  and  3  show  the  fits  to  the  elastic  scattering  cross- 
section  and  polarization  data  at  200  MeV.  The  data  are  from  Meyer  et 
al.  (Me81) .  The  experimental  data  were  reduced  to  60°,  since  there  is 
a  large  disagreement  in  the  fit  according  to  the  two  models  for  angles 
beyond  that.  The  fit  is  comparable  or  equivalent  to  that  of  Meyer  et 
al.  (Me81) .  The  reaction  cross-section  resulting  from  the  calculation 
of  the  standard  model  equals  292.7  milli  barns  (mb),  while  that  of  the 
Dirac  model  equals  236.5  mb. 

24 

VI. 2  Mg 

Figures  4  and  5  show  the  fits  to  the  elastic  scattering  cross- 
section  and  polarization  experimental  data  at  155  MeV.  The  data  are 


. 
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Figure  2  Differential  cross-section  for  P  +  12C  elastic 

scattering  at  Ep  =  200  MeV.  The  solid  curve  is 
the  DIRAC  model  and  the  dashed  curve  is  the 
STANDARD  model. 
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Figure  3  Polarization  for  P  +  12C  elastic  scattering  at 

Ep  =  200  MeV.  The  solid  curve  is  the  DIRAC  model 
and  the  dashed  curve  is  the  STANDARD  model. 
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Figure  4  Differential  cross-section  for  P  +  24Mg  elastic 
scattering  at  Ep  =  155  MeV.  The  solid  curve  is 
the  DIRAC  model  and  the  dashed  curve  is  the 
STANDARD  model. 
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Polarization  for  P  +  2LfMg  elastic  scattering  at 
Ep  =  155  MeV.  The  solid  curve  is  the  DIRAC  model 
and  the  dashed  curve  is  the  STANDARD  model. 


Figure  5 
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from  Willis  et  al.  (Wi68)  .  The  fits  according  to  the  two  models  are  nearly 
equivalent.  The  Dirac  model  is  slightly  in  better  agreement  with  the  ex¬ 
periment  for  the  elastic  polarization.  The  reaction  cross-section  resulting 
from  the  calculations  of  the  standard  model  equals  402.5  mb  while  that  of 
the  Dirac  model  equals  349.2  mb. 

VI. 3  4°Ca 

For  calcium  we  have  analysed  the  data  from  the  four  energies, 

155  MeV,  181  MeV,  300  MeV  and  500  MeV.  We  show  the  fits  to  the  elastic 
scattering  cross-section  and  polarization  at  181  MeV,  in  figures  6  and 
7.  The  data  are  from  Ingemarsson  et  al .  (In71)  .  As  the  figures  show  the 
agreement  between  the  two  models  in  fitting  the  experimental  data  is 
very  good.  They  are  approximately  equivalent.  The  reaction  cross- 
section  resulting  from  the  calculations  of  the  standard  model  equals 
495.5  mb  while  that  of  the  Dirac  model  equals  594.5  mb  and  the  experi¬ 
mental  value  (Jo61)  equals  524±14  mb. 

VI. 4  Potential  Equivalence 

It  is  known  that  the  non-relativistic  optical  potentials,  i.e.  the 
potentials  resulting  from  the  standard  model  are  energy  dependent.  The 
real  central  potential  remains  attractive  as  the  energy  increases,  then 
turn  repulsive  at  some  energy.  While  the  imaginary  part  remains  attrac¬ 
tive.  The  energy  dependence  comes  out  explicitly  in  the  Dirac  model 
(see  equation  (IV. 20)).  With  increasing  energy,  the  real  central  poten¬ 
tial  changes  slowly  from  attractive  to  repulsive  in  the  interior,  re¬ 
sulting  in  the  "Wine  Bottle  Bottom"  shape,  in  the  energy  range  100-200 

40 

MeV.  Figure  8  shows  the  real  central  potential  for  the  case  of  P  +  Ca 


at  181  MeV,  as  it  is  predicted  by  both  models,  illustrating  the  above 
mentioned  properties  of  this  potential.  Figure  9  shows  the  prediction 
of  the  imaginary  central  potential  for  the  same  case  mentioned  above. 
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Figure  6  Differential  cross-section  for  P  +  4 0 Ca  elastic 
scattering  at  Ep  =  181  MeV.  The  solid  curve  is 
the  DIRAC  model  and  the  dashed  curve  is  the 
STANDARD  model. 
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Figure  7  Polarization  for  P  4-  40Ca  elastic  scattering 

at  Ep  =  181  MeV .  The  solid  curve  is  the  DIRAC 
model  and  the  dashed  curve  is  the  STANDARD 
model . 
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Figure  8  The  real  part  of  the  central  potential  for 
P  +  4 °Ca  at  Ep  =  181  MeV.  The  solid  curve 
is  the  DIRAC  model  and  the  dashed  curve  is 
the  STANDARD  model. 
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Figure  9  The  imaginary  part  of  the  central  potential 
for  P  +  40Ca  at  Ep  =  181  MeV.  The  solid 
curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 
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Note  that  the  Dirac  model  prediction  is  much  deeper  than  that  of  the 
Standard  model.  In  fact  the  imaginary  central  potential  gets  deeper 
as  energy  increases. 

The  spin-orbit  potentials  also  have  some  energy  dependence,  but 
it  is  much  smaller  than  that  of  the  central  potential.  This  too  comes 
about  as  a  result  of  the  energy  dependent  denominator  of  equation 
(IV. 21).  In  figures  10  and  11,  we  show  both  the  real  and  imaginary 
spin-orbit  potentials  as  predicted  by  both  models.  We  can  see  that 
the  predictions  for  the  real  spin-orbit  potential  of  both  models  have  the 
same  general  shape  with  only  a  slight  difference  in  magnitude.  For 
the  imaginary  part,  on  the  other  hand,  the  difference  in  magnitude  gets 
larger.  The  Standard  model  prediction  is  higher  than  that  of  the  Dirac 
model  by  a  factor  of  nearly  3. 

Finally  one  can  say  that  the  two  potentials  (the  Standard  and  the 
Dirac  models)  do  differ  in  shape,  particularly  near  200  MeV,  but  inspite 
of  this  fact  they  still  can  be  regarded  as  equivalent,  since  they  give 
reasonably  good  agreement  with  the  data  at  this  energy  range.  This 
accord  may  not  hold  beyond  this  energy  range  due  to  the  inability  of 
the  Standard  model  to  provide  an  equivalent  fit  to  the  experimental 


data  as  that  of  the  Dirac  model. 
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Figure  10  The  real  part  of  the  spin  orbit  potential  for 
P  +  40Ca  at  Ep  =  181  MeV.  The  solid  curve  is 
the  DIRAC  model  and  the  dashed  curve  is  the 
STANDARD  model. 
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Figure  11  The  imaginary  part  of  the  spin  orbit  potential 
for  P  +  40Ca  at  E  =  181  MeV.  The  solid  curve 
is  the  DIRAC  model  and  the  dashed  curve  is  the 
STANDARD  model. 


VII.  INELASTIC  SCATTERING 


In  the  preceding  chapter  we  have  shown  comparisons  between  the  Stan¬ 
dard  model  and  the  Dirac  model  for  some  elastic  scattering  data.  In  this 
chapter  we  shall  carry  out  a  discussion  about  the  extended  optical  model 
potential  and  its  applications  in  inelastic  scattering.  In  particular ,  we 
are  interested  in  finding  out  whether  the  near  equivalence  of  the  T  Standard' 
and  'Dirac'  models  for  elastic  scattering  does  extend  to  the  case  of  inelas¬ 
tic  scattering.  In  other  words ,  we  wish  to  investigate  the  impact,  on 
inelastic  scattering,  of  the  differences  in  potentials  associated  with 
the  two  models,  particularly  in  the  region  near  200  MeV. 

We  use  the  extended  Standard  model  in  a  non -relativistic  distorted 
wave  Born  ap p r ox ima tion(DWBA)  calculation  for  proton  inelastic  scat¬ 
tering.  We  have  used  a DWBA  code  (Sh68)  for  this  purpose,  which  is  based 
on  the  formalism  given  in  chapter  III.  The  only  adjustable  parameter  in 
this  calculation  is  the  deformation  parameter  (3) .  We  have  used  equal 
deformation  length  (3R)  in  our  calculations,  since  it  gives  better  pre¬ 
diction  for  the  polarization  than  using  equal  deformation  parameters 
(see  below) . 

We  have  modified  the  (DWBA)  code  mentioned  above  to  carry  out  simi¬ 
lar  calculation  for  inelastic  scattering  based  on  the  Dirac  model. 

Three  more  subroutines  were  added  to  the  code  to  calculate  the  relati¬ 
vistic  kinematics,  the  Dirac  equation  based  "equivalent"  potentials  and 
their  associated  form  factors.  The  formalisms  of  these  calculations 
are  given  in  chapters  IV  and  V.  We  have  also  used  equal  deformation 
lengths  for  all  parts  of  the  potentials  for  the  same  reason  mentioned 
above.  For  the  purpose  of  comparisons  between  the  two  models  we  have 
used  the  same  deformation  length  in  both  cases. 
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Figures  12  to  17  show  comparisons  of  the  inelastic  scattering 

40 

form  factors  for  the  excitation  of  the  3  state  at  3.73  MeV  in  Ca  at 
181  MeV.  It  is  clear  that  for  the  real  central  form  factor  (figure  12) 
there  is  a  difference  between  the  predictions  of  the  two  models  in  the 
interior,  reflecting  the  Wine  Bottle  Bottom  shape  in  the  Dirac  model. 
Figure  13  shows  that  the  imaginary  central  form  factor  is  much  stronger 
in  the  Dirac  model  but  has  the  same  shape  except  at  very  short  distance. 
This  feature  could  be  explained  by  relating  it  to  the  fact  that  the 
imaginary  central  potential  in  the  Dirac  model  is  much  deeper  than  that 
of  the  Standard  model. 

The  spin-orbit  potential  was  deformed  and  included  in  the  calcula¬ 
tion.  In  the  present  analysis  we  use  the  full  Thomas  form  (Sh68a)  for 
the  deformed  spin-orbit  potential.  Spin-orbit  1  refers  to  the  spin- 
orbit  term  in  equation  (III. 27)  and  its  counterpart  in  equation  (V.29), 
while  spin-orbit  2  refers  to  the  spin -orbit  term  in  equation  (III. 2 8) 
and  its  counterpart  in  equation  (V.30).  We  note  from  figures  (14)  to 
(17)  that  the  spin-orbit  form  factors  in  the  two  models  are  similar  in 
shape  but  vary  in  magnitude  near  peaks  and  minima. 

We  have  found  that  the  effect  of  including  spin-orbit  deformation 
in  the  transition  potential  on  the  calculated  cross-section  is  impor¬ 
tant.  This  is  contrary  to  what  Satchler  (Sa83)  had  expected;  namely 
that  the  effect  will  not  be  appreciable  in  the  energy  range  near  200 
MeV.  Figures  18  and  19  show  an  illustration  for  the  above  mentioned 

effect  calculated  using  both  the  Standard  and  the  Dirac  models  for  the 

40 

case  of  inelastic  scattering  on  Ca  at  181  MeV.  We  can  see  that  the 
differences  are  appreciable,  they  are  present  both  in  the  magnitude  of 
the  peak  in  the  cross-section  and  in  the  angular  distribution  itself. 
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Figure  12  The  real  part  of  the  central  radial  form  factor 
for  P  +  40Ca  at  Ep  =  181  MeV.  The  solid  curve 
is  the  DIRAC  model  and  the  dashed  curve  is  the 
STANDARD  model. 
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Figure  13  The  imaginary  part  of  the  central  radial  form 

factor  for  P  +  ^Ca  at  Ep  =  181  MeV.  The  solid 
curve  is  the  DIRAC  model  and  the  dashed  curve  is 
the  STANDARD  model. 
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Figure  14  The  real  part  of  spin  orbit  1  form  factor 
for  P  +  l+0Ca  at  Ep  =  181  MeV.  The  solid 
curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 
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Figure  15  The  imaginary  part  of  spin  orbit  1  form  factor  for 
P  +  ^Ca  at  Ep  =  181  MeV.  The  solid  curve  is  the 
DIRAC  model  and  the  dashed  curve  is  the  STANDARD 
model. 
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Figure  16  The  real  part  of  spin  orbit  2  form  factor 
for  P  +  40Ca  at  Ep  =  181  MeV .  The  solid 
curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 
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Figure  17  The  imaginary  part  of  spin  orbit  2  form 
factor  for  P  +  40Ca  at  Ep  =  181  MeV. 

The  solid  curve  is  the  DIRAC  model  and 
the  dashed  curve  is  the  STANDARD  model. 
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The  effect  of  the  deformed  spin  orbit  term  on  the 
calculated  cross-section.  The  calculations  are 
the  standard  model  predictions  for  40Ca  (3“  state, 
Ex  =  3.73  MeV)  at  Ep  =  181  MeV. 
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Figure  19  The  effect  of  the  deformed  spin  orbit  term  on  the 

calculated  cross-section.  The  calculations  are 
the  DIRAC  model  predictions  for  t+0Ca  (3“  state, 

Ex  =  3.73  MeV)  at  Ep  =  181  MeV. 
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The  prediction  of  the  two  models  for  the  cross-section  drops  off  more 
rapidly  and  oscillates  more  sharply  as  the  angle  increases.  Besides 
these  observations,  it  is  known  that  the  spin-orbit  deformation  is 
important  and  crucial  in  the  calculation  of  inelastic  polarization. 

Below  we  will  discuss  some  other  effects  on  the  calculations  namely 
the  effects  of  using  equal  deformation  length,  the  effects  of  the  kine¬ 
matics  and  the  effects  of  Coulomb  deformations. 

VII. 1  The  Effect  of  Using  Equal  Deformation  Length. 

We  have  done  the  calculations  of  the  Standard  model  with  equal 
deformation  parameter  3,  which  means  all  parts  of  the  form  factors  get 
multiplied  by  the  same  number  3,  then  we  repeated  the  whole  calcula¬ 
tions  with  equal  deformation  length  3Rj  which  means  that  all  parts  of 
the  form  factors  get  multiplied  by  3Rq  where  Rq  is  the  radius  of  the 
real  part  of  the  central  potential.  We  have  found  that  the  prediction 
of  the  Standard  model  for  the  proton  inelastic  polarization  is  in  better 
agreement  with  experimental  data  when  equal  deformation  length  3R  is 
used.  So  we  have  considered  the  calculations  of  both  models  with  equal 
deformation  length  3R  in  our  analysis,  and  we  kept  3R  equal  for  both 
models  for  purposes  of  comparison. 

VII. 2  The  Effect  of  Kinematics. 

We  have  studied  the  role  of  kinematics  in  the  close  accord  between 
the  predictions  of  the  two  models.  We  have  carried  out  the  entire 
Standard  model  calculations  with  non-relativistic  and  with  relativistic 
kinematics.  We  have  compared  the  results  and  we  found  that  there  are 
minor  differences  and  one  can  say  that  the  two  types  of  calculations 
are  basically  the  same  at  the  energy  range  of  200  MeV ,  aside  from  a 


. 


. 


64 


Figure  20  Differential  cross-section  for  P  +  12C  inelastic 
scattering  (2+,  Ex  =  4.43  MeV)  at  Ep  =  200  MeV . 
The  solid  curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 


ANALYSING  POWER 


200  MEV 

-  STANDARD 

-  DIRAC 


65 


12  c 


o 


ANGLE  DEGREES) 


Figure  21  Analysing  power  for  P  4-  C  inelastic  scattering 
(2“  state,  Ex  =  4.43  MeV)  at  Ep  =  200  MeV . 

The  solid  curve  is  the  DIRAC  model  and  the 
dashed  curve  is  the  STANDARD  model. 
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Figure  22  Differential  cross-section  for  P  +  24Mg  inelastic 

scattering  (2+  state,  =  1.37  MeV)  at  Ep  =155  MeV. 
The  solid  curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 
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Figure  23  Polarization  for  P  +  24Mg  inelastic  scattering 
(2+  state,  Ex  =  1.37  MeV)  at  Ep  =  155  MeV .  The 
solid  curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 
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slight  change  in  normalization. 

VII. 3  The  Effect  of  Coulomb  Deformation. 

We  have  also  studied  the  effects  of  including  the  deformation  of 

the  Coulomb  potential,  and  we  have  found  that  the  effects  are  very 

208 

small  and  do  not  affect  the  calculations  except  for  Pb.  We  have, 
therefore,  left  these  effects  out  of  the  calculations  for  both  models. 

Below  we  will  show  comparisons  between  the  predictions  of  the 
Dirac  model  and  the  Standard  model,  as  well  as  compare  these  to  the 
experimental  results  whenever  available. 

VII. 4  12C(p,p’)12c''  2+  State  (Ex  =  4.43  MeV)  at  Ep  =  200  MeV. 

Figures  20  and  21  show  the  prediction  of  the  two  models  for  both 
inelastic  cross-section  and  analysing  power.  The  agreement  between  the 
two  types  of  calculations  is  satisfactory.  The  agreement  with  the  ex¬ 
perimental  data  is  not  that  good  but  it  is  comparable  with  other  ana¬ 
lyses  (Me81) .  The  Dirac  model  is  in  slightly  better  agreement  with  the 
data  especially  at  forward  angles.  The  data  are  from  (Me81) .  We  have 
normalized  the  Dirac  calculations  to  the  experimental  data,  then  we 
extracted  3=0.60,  and  0.439  for  the  Dirac  model  and  the  Standard  model, 
respectively,  and  3R  =  1.30  F  for  both  models.  If  we  normalize  the 
Standard  model  calculations  to  the  data  we  get  3=0*56  and  0.41  for  the 
Dirac  model  and  the  Standard  model,  respectively  with  3R  =  1*21  F  for 
both  models. 

VII.  5  2Slg(p  ,p f )  2ZfMg  2+  State  (Ex  =1.37  MeV)  at  Ep  =  155  MeV. 

Figures  22  and  23  show  comparisons  between  the  predictions  of  the 
Standard  model  and  the  Dirac  model  for  inelastic  scattering  cross- 
section  and  polarization,  and  the  experimental  data.  The  data  are  from 
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(Wi68) .  The  agreement  between  the  prediction  of  the  two  models  is 
satisfactory.  They  do  predict  the  general  shape  of  the  experimental 
data.  We  have  used  equal  deformation  length  in  both  models,  (3R  = 1.364 
fm)  ,  where  3=0.34  for  the  Standard  model  while  3  =0.482  for  the  Dirac 
model.  Sherif  (Sh68) ,  in  his  analysis  of  the  same  case  using  the  Stan¬ 
dard  model  obtained  3  =0.353  and  3R  =  1.30  fm. 

32  3?  *  + 

VII. 6  S(p,pT)  S  2  State  (Ex  =2.24  MeV)  at  Ep  =  155  MeV. 

Figures  24  and  25  show  comparisons  between  the  predictions  of  the 
Standard  and  the  Dirac  models  for  the  inelastic  cross-section  and  pola¬ 
rization.  We  also  show  the  experimental  data  of  Willis  et  al.  (Wi68) 
for  comparison.  The  agreement  between  the  two  types  of  calculations  is 
very  good.  Where  the  agreement  between  the  calculations  and  the  experi¬ 
mental  data  is  quite  good  for  the  inelastic  cross-section  over  the  whole 
angular  range  (i.e.  up  to  0=60°)  while  it  is  good  for  the  polarization 
only  up  to  0=40°.  We  get  a  deformation  parameter  3=0.23  for  the 
Standard  model  and  3  =  0.31  for  the  Dirac  model  and  a  3R  =  1.033  F  for 
both  models.  From  previous  analyses  for  the  same  case,  Sherif  (Sh68) 
obtained  3  =  0.24  and  3R  =  1.15  F  while  Haybron  (Ha66)  got  a  fit  to  an¬ 
other  set  of  data  using  collective  model  DWBA  without  spin-orbit  defor¬ 
mation  and  obtained  3=0.29  and  3R  =  1.12  F. 

40  40  * 

VII.  7  Ca(p ,p 1 )  Ca  3  State  (Ex  =  3.73  MeV)  at  Ep  =  181  MeV. 

Figures  26  and  27  show  the  predictions  of  the  two  models  for  both 
inelastic  cross-section  and  polarization  comparing  them  with  the  experi¬ 
mental  data.  The  data  are  from  (In71) .  The  agreement  between  the  pre¬ 
dictions  of  the  two  models  is  very  good.  The  results  of  the  present 
work  are  superior  to  that  of  Satchler  (Sa83)  in  fitting  the  experimental 
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Figure  24  Differential  cross-section  for  P  +  32S  inelastic 
scattering  (2+  state,  Ex  =  2.24  MeV)  at  Ep  =  155 
MeV.  The  solid  curve  is  the  DIRAC  model  and  the 
dashed  curve  is  the  STANDARD  model. 
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Figure  25  Polarization  for  P  +  32S  inelastic  scattering 

(2+  state,  Ex  =  2.24  MeV)  at  Ep  =  155  MeV.  The 
solid  curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 


' 


72 


ANGLE  (DEGREES) 


Figure  26  Differential  cross-section  for  P  +  40Ca  inelastic 

scattering  at  Ep  =  181  MeV.  The  solid  curve  is 
the  DIRAC  model  and  the  dashed  curve  is  the  STANDARD 
model. 
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Polarization  for  P  +  40Ca  inelastic  scattering 
(3“  state,  Ex  =  3.73  MeV)  at  Ep  =  181  MeV.  The 
solid  curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 


Figure  27 
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data,  since  he  neglected  the  effects  of  the  deformed  spin-orbit  poten¬ 
tial  which  were  shown  to  be  important  for  this  energy  range.  We  get  a 
deformation  parameter  3=0.23  for  the  Standard  model,  3  =0.285  for  the 
Dirac  model  and  hence  3R  = 1.014  F  for  both  models.  Satchler  (Sa83) 
used  equal  deformation  parameter  3=0.35  in  his  analysis  and  3R  =  1.464 
and  1.398  F  for  the  W-S  model  and  the  WBB  model  respectively.  If  we 
carry  out  our  calculations  without  including  the  effect  of  the  deformed 
spin-orbit  interaction  we  get  3=0.27  and  3  =0.336  for  the  Standard 
model  and  the  Dirac  model  respectively  and  3R  =  1.19  F  for  both  of  them. 
The  ratio  of  the  peak  cross-section  of  the  Standard  model  to  that  of 
the  Dirac  model  with  and  without  including  the  spin-orbit  interaction 
in  the  calculation  is  1.093  and  1.168  respectively  if  we  use  equal  de¬ 
formation  length.  If  we  use  equal  deformation  parameter  3=0.27  for 
the  calculations  without  including  the  spin-orbit  the  ratio  will  be 
1.81,  where  according  to  Satchler  (Sa83)  the  ratio  will  be  2  if  we  use 
equal  deformation  parameter  (3),  and  1.5  if  we  use  equal  deformation 
length  (3R) . 

40  4f)  *  + 

VII. 8  Ca(p ,p ’ )  Ca  2  State  (Ex  = 18  MeV)  at  Ep  = 181  MeV  . 

Figures  28  and  29  show  the  predictions  of  both  models  for  inelas¬ 
tic  scattering  and  polarization.  We  don’t  have  experimental  data  for 
this  case  to  compare  with  the  theoretical  calculations,  but  we  note 
that  the  two  models  are  in  good  agreement.  We  have  used  3=0.23  for 
the  Standard  model  and  3  =0.285  for  the  Dirac  model  and  3R  =  1.014  F  for 
both  models.  Satchler  (Sa83)  used  3  =  0.254  for  both  W-S  and  WBB  models 
and  3R  =  1»094  F  for  the  W-S  model  and  0.981  for  the  WBB  model.  The 
peak  cross-section  for  the  Standard  model  is  15.57  mb/sr  while  it  is 
12.75  mb/sr  for  the  Dirac  model,  Satchler  (Sa83)  obtained  a  value  of 
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Figure  28  Differential  cross-section  for  P  +  40Ca  inelastic 
scattering  (2+  state,  Ex  =  18  MeV)  at  Ep  =  181 
MeV.  The  solid  curve  is  the  DIRAC  model  and  the 
dashed  curve  is  the  STANDARD  model. 
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Figure  29  Polarization  for  P  +  Lf0Ca  inelastic  scattering 

(2+  state,  Ex  =  18  MeV)  at  Ep  =  181  MeV.  The 
solid  curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 
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14.9  mb/sr  for  the  W-S  model  and  7.9  mb/sr  for  the  WBB  model. 

40  40  *  - 

VII. 9  Ca(p,pf)  Ca  3  State  (3.73  MeV)  at  Ep  =  300  MeV. 

Figures  30  and  31  show  the  predictions  of  both  models  for  the  in¬ 
elastic  cross-section  and  polarization.  The  agreement  between  the  two 
models  is  good  in  cross-section,  while  large  differences  appear  in  po¬ 
larization,  There  is  no  experimental  data  for  this  case  to  test  the 
two  models  T predictions .  It  should  be  mentioned  here  that  the  two  models 
do  not  predict  the  same  elastic  scattering,  the  Dirac  model  is  superior 
to  the  Standard  model.  We  believe  that  the  differences  in  the  inelastic 
scattering  relations  are  due  to  those  of  the  elastic  scattering. 

40  40  *  - 

VII. 10  Ca(p ,p ’ )  Ca  3  State  (3.73  MeV)  at  500  MeV. 

Figures  32  and  33  show  the  predictions  of  the  two  models  for  both 
inelastic  cross-section  and  polarization.  There  are  no  published  ex¬ 
perimental  data  for  this  specific  case  but  some  preliminary  experimental 
data  have  been  made  available  to  us  (Se83) .  At  this  energy  the  dis¬ 
agreement  between  the  predictions  of  the  two  models  is  getting  larger, 
particularly  for  the  polarization.  It  should  be  mentioned  here  that 
the  two  models  do  not  predict  the  same  elastic  scattering  at  this  energy, 
the  Dirac  model  is  superior  to  the  Standard  model.  We  believe  that  the 
differences  in  the  predictions  of  the  inelastic  scattering  observables 
are  a  magnification  of  those  of  the  elastic  scattering. 

208  OQQ  ^  _ 

VII. 11  Pb(p,p?)  Pn  3  State  (2.62  MeV)  at  Ep  = 200  MeV. 

Figures  34  and  35  show  the  predictions  of  the  two  models  for  in¬ 
elastic  scattering  cross-section  and  polarization.  The  experimental 
data  are  from  (Be 81) .  The  two  models  are  in  good  agreement  in  the  cross- 
section  predictions,  while  they  are  in  less  agreement  in  the  polarization 
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Figure  30  Differential  cross-section  for  P  +  ^Ca 

inelastic  scattering  (3“  state,  Ex  =  3.73  MeV) 
at  Ep  =  300  MeV.  The  solid  curve  is  the  DIRAC 
model  and  the  dashed  curve  is  the  STANDARD  model. 
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Polarization  for  P  +  40Ca  inelastic  scattering 
(3“  state,  Ex  =  3.73  MeV)  at  Ep  =  300  MeV. 

The  solid  curve  is  the  DIRAC  model  and  the 
dashed  curve  is  the  STANDARD  model. 
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Figure  32  Differential  cross-section  for  P  4-  40Ca  inelastic 

scattering  (3“  state,  Ex  =  3.73  MeV)  at  Ep  =  500 
MeV.  The  solid  curve  is  the  DIRAC  model  and  the 
dashed  curve  is  the  STANDARD  model. 
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Figure  33  Polarization  for  P  +  40Ca  inelastic  scattering 

(3“  state,  Ex  -  3.73  MeV)  at  Ep  =  500  MeV .  The 
solid  curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 
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Figure  34  Differential  cross-section  for  P  +  288Pb 

inelastic  scattering  (2+  state,  Ex  =  2.62 
MeV)  at  Ep  =  200  MeV.  The  solid  curve  is 
the  DIRAC  model  and  the  dashed  curve  is 
the  STANDARD  model. 
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Polarization  for  P  +  208Pb  inelastic  scattering 
(2+  state,  Ex  =  2.62  MeV)  at  Ep  =  200  MeV.  The 
solid  curve  is  the  DIRAC  model  and  the  dashed 
curve  is  the  STANDARD  model. 


Figure  35 
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prediction.  We  used  3=0.095  for  the  Standard  model  and  3=0.1132  for 
the  Dirac  model  and  3R  =  0.745  F  for  both  models.  Satchler  (Sa83)  used 
3=0.116  in  his  analysis  for  both  models,  and  3R  =  0 . 8405  fm  for  the 
W-S  model  and  3R= 0.8027  fm  for  the  WEB  model.  He  did  not  include  the 
effect  due  to  the  deformed  spin-orbit  potential  which  we  found  is  im¬ 
portant  and  necessary.  The  prediction  for  peak  cross-section  is  (23.8 
mb/sr)  for  the  W-S  model,  while  that  of  the  WBB  is  16.70  mb/sr  which 
means  that  the  ratio  between  the  predictions  of  the  two  models  is  1.425. 
In  our  analysis  the  Standard  model  predicts  a  peak  cross-section  of 
23.4  mb/sr  while  the  Dirac  model  predicts  16.84  mb/sr,  this  means  that 
the  ratio  between  the  predictions  of  the  two  models  is  1.39.  If  we  do 
the  calculations  without  including  the  deformed  spin-orbit  interaction 
the  Dirac  model  predicts  a  peak  cross-section  of  14.675  mb/sr  using 
the  same  3k  as  before,  but  if  we  fit  the  calculations  to  the  experimen¬ 
tal  data  we  get  3=0.143  and  3R  =  0.94  F.  The  shape  of  the  angular  dis¬ 
tribution  will  be  different,  the  effect  is  illustrated  in  Figure  36. 
From  the  above  discussion,  it  is  clear  that  the  effect  of  using  the 
deformed  spin  dependent  interaction  is  important  especially  in  predic¬ 
ting  the  shape  of  the  angular  distribution  (see  figure  36).  It  has 

less  effect  in  the  prediction  of  the  magnitude  of  the  peak  cross-section, 

40 

than  the  case  of  3  state  in  Ca  at  Ep  =  181  MeV. 

The  effects  of  Coulomb  deformations  for  this  case  have  been  inves¬ 
tigated.  We  have  done  the  calculations  of  the  Standard  model  again, 
but  including  the  deformed  Coulomb  potential.  We  have  found  that  the 
Standard  model  predicts  a  peak  cross-section  of  20.31  mb/sr,  which 
means  that  the  inclusion  of  the  deformed  Coulomb  potential  does  reduce 
the  peak  cross-section  by  about  15%.  The  effects  are  not  large  and  they 
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Figure  36  The  effect  of  the  deformed  spin  orbit  term  on 

the  calculated  cross-section.  The  calculations 
are  the  Dirac  model  predictions  for  288Pb  (3” 
state,  Ex  =2.61  MeV)  at  Ep  = 200  MeV. 
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The  effect  of  the  deformed  Coulomb  potential  on 
the  calculated  cross-section.  The  calculations 
are  the  Standard  model  predictions  for  288Pb  (3” 
state.  Ex  =2.61  MeV)  at  Ep  =  200  MeV. 


Figure  37 
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do  not  change  the  shape  of  the  distribution.  The  effect  is  illustrated 
in  Figure  37. 

Satchler  (Sa83)  calculates  the  values  of  3  for  some  low-lying 
collective  states  from  electromagnetic  transitions,  but  his  values  for 
3  were  higher  than  those  of  the  present  work.  We  believe  that  one  of 
the  reasons  for  this  is  that  Satchler  (Sa83)  did  not  include  the  ef¬ 
fects  of  the  deformed  spin-orbit  interaction  in  his  analysis  which  we 
found  would  require  higher  3  to  fit  the  experimental  data.  Osterfeld 

et  al.  (Os79) ,  in  a  recent  investigation  of  the  hadronic  transition 

208  — 

operator  for  inelastic  proton  scattering  from  Pb  (3  state.  Ex  =2.61 
MeV) ,  have  found  that  the  value  of  3R  calculated  by  that  operator  de¬ 
viates  about  24%  from  that  calculated  for  the  electromagnetic  transir- 
tion  operator  at  Ep  = 40  MeV,  while  the  deviation  reaches  about  50%  for 
the  same  case  at  Ep  = 156  MeV.  They  also  mentioned  that  the  value  of 
3R  decreases  as  the  incident  projectile  energy  increases.  As  we  can 
see  from  table  III,  Osterfeld  et  al .  (Os79)  obtained  these  values  for 
3R  =  0 . 899 ,  0.735,  0.638,  0.474  F  at  Ep  =24.55,  40.0,  54.0,  61.25,  156.0 
MeV,  respectively.  Satchler  (Sa83)  has  obtained  R= 0.8405,  0.8027  F 
for  the  W-S  model  and  the  WBB  model  respectively  at  Ep=200  MeV  for 
the  same  case.  In  the  present  work  we  have  obtained  a  3R  = 0.745  F  for 
both  models.  Gruhn  et  al.  (Gr72)  also  noticed  the  same  feature  of  de¬ 
creasing  3R  as  projectile  energy  increases  in  their  analysis  of  inelas- 

40 

tic  proton  scattering  from  Ca  at  different  projectile  energies.  They 

have  obtained  3R  =  1*44,  1.40,  1.32,  1.32  F  for  proton  inelastic  scatter- 
40 

ing  from  Ca  (3  state;  Ex  =  3.73  MeV)  at  Ep=17,  25,  40,  55  MeV  res¬ 
pectively.  Sherif  (Sh68)  obtained  3R = 1.187  F  at  Ep  = 155  MeV  for  the 
same  case.  Satchler  (Sa83)  got  3R  = 1.464,  1.398  F  for  the  W-S  and  the 


TABLE  III 


THE  VALUES  OF  $R(F)  FOR  INELASTIC  PROTON 
208  — 

SCATTERING  FROM  Pb  (3  STATE;  Ex  =  2.61  MeV) 


Ep  (MeV) 

Osterfeld  et  al. 
(Os72) 

Satchler 

(Sa83) 

Present  Work 

24.55 

0.899 

— 

— 

40.0 

0.725 

— 

— 

54.0 

0.732 

— 

— 

61.25 

0.638 

— 

— 

156.0 

0.474 

— 

200.0 

... 

0.8405  (W-S) 

0.745 

0.8027  (WBB) 

TABLE  IV 


THE  VALUES  OF  $R(F)  FOR  INELASTIC  PROTON 
40 

SCATTERING  FROM  Ca  (3  STATE;  Ex  =  3.  73  MeV)  . 


Ep  (MeV) 

Gruhn  et  al. 
(Gr  72) 

Sherif 
(Sh  68) 

Satchler 

(Sa83) 

Present  work 

17.0 

1.44 

— 

— - 

— 

23.0 

1.40 

— 

— 

40.0 

1.32 

— 

— 

55.0 

1.32 

— 

— 

— 

155.0 

— 

1.187 

— 

1.058 

181.0 

1.464  (W-S) 

1.014 

1.398  (WBB) 

WBB  models  at  Ep  = 181  MeV  respectively.  In  the  present  work  we  have 
obtained  3R  =  1.014  F  for  both  models  at  Ep  =  181  MeV. 

Finally,  we  would  like  to  say  that  the  values  of  3R  obtained  in 
the  present  work  do  agree  with  the  general  behaviour  of  (SR>  mentioned 
in  the  above  discussion,  more  than  those  of  Satchler  (Sa83) .  Inspite 
of  that,  the  magnitude  of  (3  remains  an  open  question. 


VIII.  CONCLUSION 


The  work  described  in  this  thesis  was  done  to  study  both  conven¬ 
tional  and  unconventional  optical  potential  models  for  proton  inelastic 
scattering.  We  have  carried  out  comparisons,  over  a  wide  range  of  bom¬ 
barding  energies,  between  the  predictions  of  two  models,  namely  the 
Standard  model  based  on  Schrodinger  equation  and  the  Dirac  model  based 
on  the  reduced  Dirac  equation,  and  experimental  data  whenever  available 
at  the  same  or  nearby  energy. 

We  have  found  in  the  present  work  that  the  deformed  spin-dependent 
interactions  in  addition  to  being  crucial  in  the  calculation  of  polari¬ 
zation,  are  relatively  important  in  the  case  of  cross-section,  in  the 
energy  region  of  155-500  MeV  .  We  have  carried  out  calculations  with  and 
without  including  the  deformed  spin -dependent  interactions  in  this  energy 
range  and  we  have  found  that  they  are  important  for  predicting  both  the  right 
magnitude  for  the  peak  cross-section  and  the  general  shape  of  the  ex¬ 
perimental  data.  We  also  concluded  that  the  kinematics  do  not  crucial¬ 
ly  affect  the  results  at  this  energy  aside  from  a  slight  change  in 
normalization . 

We  have  used  equal  deformation  lengths  for  the  different  parts  of 
the  potentials  throughout  this  work,  and  we  have  found  that  this  is  an 
important  and  a  necessary  ingredient  for  the  calculations,  because  the 
results  obtained  for  the  inelastic  polarization  are  in  a  better  agree¬ 
ment  with  the  data  than  those  using  equal  deformation  parameters  (3) • 

We  would  like  to  emphasize  here  that  the  change  in  the  shape  of 
the  potential  does  not  crucially  affect  the  calculation  if  all  parts  of 
the  potential  including  the  spin-dependent  part  are  deformed  and  included 
in  a  full  consistent  manner,  with  a  necessary  constraint  that  all  of 
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these  parts  have  the  same  deformation  length.  This  is  contrary  to 
Satchler’s  calculation  (Sa83)  ,  where  he  found  that  the  ratio  between 
the  peak  cross-section  of  Standard  to  Wine  Bottle  Bottom  models  pre¬ 
dictions  is  about  1.5  when  equal  deformation  lengths  are  used  while  it 
is  about  2  when  equal  deformation  parameters  are  used.  Our  results 
show  that  the  two  models  are  approximately  equivalent  in  their  predic¬ 
tion  of  the  inelastic  cross-section  at  this  energy  in  spite  of  the  fact 
that  the  shapes  of  the  real  central  potentials  are  quite  different  in¬ 
side  the  nuclear  volume.  We  have  also  compared  the  predictions  of  the 
two  models  for  the  inelastic  polarization.  Although  they  differ  some¬ 
what  they  still  have  the  same  general  shape. 

Such  accord,  however,  may  not  hold  as  the  projectile  energy  is 
increased.  This  would  be  mainly  due  to  the  inability  of  the  Standard 
model  to  provide  as  good  a  representation  of  the  elastic  data  as  the 
Dirac  model.  In  other  words,  the  two  potentials  are  not  expected  to  be 
nearly  equivalent  except  perhaps  near  1  GeV  bombarding  energy  where 
both  models  yield  completely  repulsive  real  central  potentials. 

The  quality  of  the  fits  to  the  angular  distribution  predicted  by 
the  Dirac  model  is  either  comparable  or  slightly  in  a  better  agreement 
with  the  experimental  data  to  those  of  the  Standard  model.  The  agree¬ 
ment  between  the  two  calculations  for  the  inelastic  scattering  observ¬ 
ables  depends  on  that  of  the  elastic  scattering.  We  believe  that  the 
differences  between  the  two  calculations  are  somehow  a  magnification  of 
their  differences  in  fitting  the  elastic  data,  which  means  that  the 
two  potentials  are  not  entirely  equivalent.  Since  the  Dirac  model  fits 
elastic  scattering  better  than  the  Standard  model  or  at  least  the  same, 
it  is  expected  to  do  the  same  in  the  inelastic  scattering  calculation 
and  this  is  what  we  conclude  from  the  present  study. 


' 
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